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ABSTRACT 


A  new  ray  tracing  computer  program  which  optimizes  the  design  of  linearly 
dispersive,  two-part  magnetoelastic  delay  lines  operating  at  microwave 
frequencies  has  been  developed.  The  simulation  shows  that  delay  line 
bandwidth  and  dynamic  range  can  be  separately  optimized  for  the  value 
of  a  wave  focusing  parameter  within  a  limited  range.  Whereas  previous 
work  had  demonstrated  that  there  are  a  potentially  infinite  number  of 
synthesized  magnetic  fields  which  produce  the  desired  linear  delay 
function,  choice  of  magnetic  field  can  now  be  made  to  obtain  a  desired 
value  of  the  wave  focusing  parameter  and  thereby  optimize  delay  line 
performance . 

The  computer  simulacion  is  correlated  with  experimental  magnetoelastic 
delay  line  results.  Previously,  linearly  dispersive  two-port  magneto¬ 
elastic  delay  lines  with  synthesized  magnetic  fields  had  been  built  in 
both  a  Strauss  type  configuration  and  with  a  novel,  "buried-antenna" 
design  whore  the  antennas  were  wires  inserted  through  separate  holes 
drilled  by  laser  in  a  YIG  rod.  The  validity  of  the  ray  tracing  is 
confirmed  by  the  experimental  data  and,  additionally,  the  simulation 
elucidates  the  causes  of  the  differences  in  bandwidth  performance 
between  the  design  types. 

T.he  critical  factors  in  optimizing  bandwidth/dynamic  range  are  shown 
to  be  the  axial  and  radial  magnetic  field  gradients  which  control  both 
the  focusing  parameter  and  the  time  evolution  of  wavelength.  A  simple 
predictive  model  for  bandwidth  is  developed  as  well  as  a  closed  form 
solution  for  the  magnetic  field  which  maximizes  bandwidth. 
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PREFACE 


The  increasing  cost  and  complexity  of  microwave  signal  generating 
and  processing  systems  continues  to  spur  interest  in  the  development  of 
components  which  operate  directly  at  microwave  frequencies.  Currently, 

t 

many  microwave  signal  processing  systems  require  elaborate  frequency 
conversion  schemes  wherein  microwave  signals  are  heterodyned  down  to 
VHF/UHF  and  are  then  processed  by  a  wide  variety  of  existing  RF  signal 
processing  devices.  The  explosive  growth  in  Surface  Acoustic  Wave  (SAW) 
devices,  for  example,  is  a  direct  result  of  the  advent  of  large  microwave 
systems  utilizing  downconvers ion  to  RF .  These  downconvers ion  methods, 
however,  can  be  extraordinarily  expensive  (as  when  extremely  coherent 
downconve rs ion  is  required)  and  lack  the  extremely  wide  bandwidths  which 
are  a  fundamental  design  objective  in  microwave  system  design.  A  case  in 
point  is  wideband  linear  FM  modulated  (CHIRP)  radar  which  suffers  from 
hardware  bandwidth  limitations  even  today. 

Magnetostatic  and  magnetoelastic  wave  devices  hold  the  promise  of 
addressing  some  of  these  issues.  Components  utilizing  the  propagation 
and  conversion  of  magnetostatic  and  acoustic  waves  offer  the  advantages 
of  direct  operation  at  microwave  frequencies,  nonreciprocal  spatial 
propagation,  relatively  long  delay  times,  and  frequency  dispersion. 
Although  many  common  devices  utilize  these  and  other  magnetostatic 
effects  (i.e.,  YIG  filters  and  oscillators,  isolators,  circulators, 
etc.),  development  of  linearly  dispersive  delay  lines  with  user 
selectable  characteristics  has  been  limited  by  difficulties  in 
synthesizing  the  requisite  internal  magnetic  fields. 


Within  the  last  10  years  a  technique  for  synthesizing  a  desired 
magnetic  field  in  a  YIG  rod  or  disc  has  become  available  [1].  The  method 
involves  the  defining  of  the  equipotential  surfaces  both  within  and 
without  the  medium.  Polepieces  which  are  equipotentials  outside  the  rod 
are  then  built  and  matched  to  the  rod;  energizing  the  pole  pieces  creates 
corresponding  internal  potential  and  magnetic  fields. 

As  a  result  of  this  field  synthesis  technique,  magnetoelastic  delay 
lines  with  linear  dispersion  can  be  made  and  optimized  to  application 
[2,3].  These  delay  lines  are  extremely  small,  and  have  attenuations  and 
bandwidths  suitable  for  most  signal  processing  applications.  This  thesis 
recapitulates  the  theory  of  operation  of  these  delay  lines,  summarizes 
the  experimental  work  performed  to  date,  and  describes  the  optimization 
technique  which  makes  delay  line  design  practicable. 


Chapter  I 


INTRODUCTION 


I .  A  Thesis  Ob-jectives 

The  goal  of  this  research  is  to  determine  whether  the  design  of  a 
ferrite  device,  the  magnetoelastic  delay  line,  can  be  improved  so  as  to 
maximize  output  bandwidth  and  dynamic  range.  The  purpose  of  the  thesis 
document  itself  is  to  evaluate  the  method,  results,  and  conclusions  of 
computer  simulation  of  that  type  of  delay  line,  and  to  describe  the 
requirements  for  successful  optimization.  The  computer  program  which 
is  the  analytical  tool  for  simulating  the  raagnetoelastic  delay  line,  the 
model  which  resulted  from  consideration  of  the  data,  and  the  conclusions 
on  optimizing  delay  line  performance  are  all  discussed  as  part  of  the 
presentat ion . 

This  introductory  chapter  is  intended  to  provide  a  perspective 
from  which  to  view  the  work  accomplished.  The  chapter  includes  a  short 
history  of  the  discovery  of  magnetostatic  and  magnetoelastic  waves  and 
their  in-jnediatc  application  in  microwave  delay  devices.  Also,  a  survey 
is  made  of  a  variety  of  magnetoelastic  devices.  Finally,  the  chapter 
concludes  with  a  somewhat  closer  look  at  the  predecessor  work  upon  which 
this  thesis  is  based. 

Chapter  II  ("Theory  of  Magnetoelastic  Waves")  and  Appendix  A  of  this 
thesis  constitute  an  abbreviated  discussion  of  the  underlying  theory  of 
magnetoelastic  waves.  Included  therein  are  a  review  of  magnetostatic  and 
magnetoelastic  wave  dispersion  relations,  wave  behaviors,  and  coupling 
effects.  Chapter  III  ("The  Magnetoelastic  Delay  Line")  and  Appendices  B 


and  C  provide  a  discussion  of  the  operation  of  the  magnetoelastic  delay 
line,  the  magnetic  field  anisotropy  of  the  YIG  rod  which  is  the  medium 
of  propagation,  and  the  field  parameters  required  for  synthesis  of  the 
desired  linear  delay  function.  Although  Chapters  II  and  III  extensively 
reference  the  appendices,  they  are  intended  to  be  a  self-contained 
introduction  sufficient  to  allow  understanding  of  the  problem  and 
subsequent  chapters. 

Chapter  IV  ("Ray  Tracing")  presents  the  theory  and  techniques  used 
in  the  ray  tracing.  Because  the  magnetic  field  in  the  ferrite  rod  is 
anisotropic,  the  propagation  of  energy  depends  upon  position  in  the 
medium  and  assumptions  more  complex  than  those  of  simple  optics  must  be 
made.  The  validity  of  the  assumptions  which  allow  ray  tracing  to  be  used 
in  an  anisotropic  medium  is  discussed  and  the  equations  which  describe 
the  time  evolution  of  the  wave  are  presented.  Also,  the  computer  program 
which  performs  the  ray  tracing  is  outlined. 

Chapter  V  ("Results  and  Conclusions")  provides  a  brief  summary  of 
the  raw  data  and  presents  conclusions  in  the  form  of  a  simple  model. 
Predictions  which  extend  the  model  to  other  substances  are  also  made,  and 
the  ultimate  goal  of  the  thesis  research,  optimizing  delay  line  design, 
is  discussed  in  the  context  of  the  model. 

Chapter  VI  ("Conclusions  and  Recommended  Future  Work")  concludes  the 
thesis  with  some  suggestions  for  future  work. 


I . B  History  of  the  Magnetoelastic  Delay  Line 


Although  spin  waves  had  been  extensively  considered  previously, 
it  was  Ernst  Schlomann  who,  in  1961,  noted  that  naturally  occurring 
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nonuniform  magnetic  fields  in  YIG  rods  and  disks  could  be  used  to  couple 
electromagnetic  energy  to  magnetostatic  spin  waves  [4].  His  observation 
was  that  long  wavelength  spin  waves  should  exist  near  surface  regions  of 
net  magnetic  dipole  moment.  They  would  then  propagate  into  the  sample 
towards  regions  of  lesser  magnetic  field  while  simultaneously  decreasing 
in  wavelength.  Eshbach  subsequently  verified  this  prediction  and  also 
noted  that  under  certain  conditions  of  magnetic  field,  the  spin  waves 
could  be  converted  into  elastic  waves  [5].  The  experimental  fact  that 
the  internal  magnetic  field  controlled  conversion  of  spin  waves  to 
elastic  waves,  and  thus  the  propagation  time,  led  to  proposals  by 
Schloraann  [6]  and  others  [7,8]  for  a  magnetoelastic  delay  line.  These 
proposals  envisioned  a  delay  line  in  which  delay  would  be  controlled  by 
tuning  the  magnitude  of  the  external  magnetic  field. 

The  delay  line  designs  of  Schlomann  and  Strauss  reflected  the 
simplicity  of  using  the  .naturally  occurring  nonuniform  internal  magnetic 
fields.  Schlomann  inserted  the  ends  of  a  YIG  rod  into  cavities. 

Strauss'  design  [7],  of  greater  utility,  mounted  a  wire  against  the 
end  face  of  a  rod  in  an  axial  magnetic  field;  the  wire  served  to 
both  transmit  and  receive  microwave  pulses.  Subsequent  designs  have 
often  been  based  upon  Strauss'  original  single  ended  conception  with 
improvements.  Both  these  designs,  however,  produced  nonlinear  dispersive 
delay  functions  over  bandwidths  typical  of  microwave  systems. 

Continued  work  in  magnetoelastic  delay  lines  has  focused  upon  two 
issues:  (a)  synthesizing  internal  magnetic  fields  which  can  provide 

linear  or  other  desired  delay  functions,  and  (b)  improving  efficiency 
and  feedthrough  isolation  between  input  and  output.  In  the  former 


case,  efforts  have  been  made  to  produce  linear  fields  by  adding  magnetic 


shunts  [9]  or  trying  unusual  crystal  shapes  [10].  These  approaches  have 
generally  proven  ineffectual.  Work  on  the  latter  has  included  improving 
the  antennas  mounted  on  rod  endfaces  with  orthogonal  thin  wire  coupling 
antennas  [2],  etched  antenna  designs  [11],  and  two-ended  polycrystalline 
YIG-YAG-YIG  constructions  [12].  Improvements  in  reproduceabil ity  and 
insertion  loss  have  been  seen. 

Borgeaud's  work  (predecessor  to  this  thesis)  in  magnetoelastic  delay 
line  manufacture  incorporated  two  extremely  effective  improvements  [3]. 
Previously,  Morgenthaler  had  developed  a  method  for  synthesizing  any 
desired  on-axis  magnetic  field  in  a  YIG  rod  or  disk  [1].  Since  the 
formulae  for  magnetic  fields  giving  linear  delay  functions  were  well 
established,  Morgenthaler  and  Platzker  were  able  to  incorporate  this 
method  in  the  construction  of  a  Strauss  type  single-ended  delay  line,  and 
they  found  much  improved  bandwidth  [2].  Their  design,  however,  still 
suffered  the  liabilities  of  single-ended  delay  line  construction. 

Borgeaud  and  Morgenthaler  then  extended  the  magnetic  synthesis 
technique  to  a  two-ended  delay  line  design,  wherein  the  antennas  were 
fine  wires  inserted  through  extremely  small  holes  drilled  in  a  YIG  rod 
by  laser.  The  synthesized  field  was  designed  to  the  requirements  of 
two-ended  operation.  Borgeaud  found  substantial  improvement  in  isolation 
and  insertion  loss  but  encountered  reduced  bandwidth  for  his  particular 
experimental  case  [3]. 

The  logical  extension  of  Borgeaud's  work  is  to  seek  greater 
bandwidth  without  sacrificing  improved  insertion  loss  and  isolation. 

There  was  no  a  priori  reason  to  expect  Borgeaud's  delay  line  to  show 
reduced  bandwidth;  the  operation  of  his  delay  line  is  identical  to  that 
of  a  single-ended  delay  line  absent  the  liability  of  a  reflection  of 


elastic  waves.  Borgeaud  had,  liowever,  selected  and  tested  only  one  of 


an  infinite  possibility  of  synthesized  magnetic  field  profiles  and  his 
selection  was  made  prior  to  the  development  of  computer  analysis  tools. 

The  purpose  of  this  thesis,  as  stated  earlier,  is  to  use  improved 
ray  tracing  techniques  to  analyze  Borgeaud's  delay  line  design.  Implicit 
in  that  is  discovering  how  to  optimize  bandwidth  without  sacrificing 
improved  isolation  and  insertion  loss,  and  explaining  the  previously 
described  results.  Development  of  a  general  model  for  delay  line 
behavior  and  consideration  of  dynamic  range  are  goals  included  within 
this  topic  as  well. 
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characterize;  it  suffices  to  say  that  electromagnetic  waves  couple  to 
MSW s  with  comparable  wavelengths  [19]. 

Acoustic  waves  in  crystal  media  are  well  understood  as  propagating 
vibrations  in  the  crystal  lattice  (Figure  2.1c)  [7].  In  fact,  magnetic 
spin  waves  and  acoustic  waves  are  inherently  coupled  together,  as  the 
magnetic  forces  produce  lattice  strain  ("magnetostriction")  and  a 
conversion  of  energy  from  magnetic  spin  wave  to  elastic  wave  [20].  Any 
propagation  of  magnetic  or  elastic  energy  in  a  f err iraagnet ic  actually 
contains  an  admixture  of  both;  thus  the  term  "magnetoelastic"  waves. 
Practically,  however,  the  proportions  of  one  or  the  other  type  of  wave 
tend  to  be  small  except  under  the  specific  conditions  which  convert  one 
wave  type  to  the  other.  This  conversion  can  be  effected  with  high 
efficiency  under  restricted  conditions  [20]. 

Spin  waves  can  thus  be  envisioned  as  a  bridge  between 
electromagnetic  and  elastic  waves.  Indeed,  spin  waves  are  an 
extraordinarily  effective  bridge  because:  (a)  they  have  wavenumbers 
which  vary  with  local  magnetic  field,  allowing  coupling  at  different 
positions  within  the  same  medium  to  modes  with  different  wavelengths; 

(b)  they  couple  to  both  EM  and  elastic  waves  with  very  high  efficiencies; 

(c)  they  have  well  defined  characteristics  which  lend  themselves  to 
theoretical  and  computer  modeling;  and  (d)  they  propagate  with  relatively 
low  attenuation.  As  analysis  will  show,  these  properties  of  spin  waves 
and  the  development  of  a  magnetic  field  synthesis  method  led  to  the 
design  of  single  ended  magnetoelastic  delay  lines  and,  ultimately,  to 
two-port  magnetoelastic  delay  lines  which  are  the  subject  of  this  work. 

The  delay  line  itself  is  a  single  crystal  rod  made  of  Yttrium  Iron 
Garnet  (Y^Fe^Oj^),  also  called  YIG.  YIG  is  f err iroagnet ic ,  having  two 
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interleaved  magnetic  sublattices  with  net  electron  spin  moments,  and  is 
optimum  for  this  application.  Its  properties  include: 

(a)  High  resistivity:  This  allows  for  penetration  of  electromagnetic 
energy  into  the  rod  and  high  efficiency  coupling  of  electromagnetic 
to  spin  energy. 

(b)  Low  saturation  magnetization:  Magnetic  domains  are  eliminated  at 
relatively  low  values  of  external  magnetic  field. 

(c)  Relatively  slow  relaxation  time:  Magnetic  and  acoustic  linewidths 
are  small  and  attenuation  from  consequent  damping  is  small. 

(d)  Easy  synthesis:  YIG  is  grown  into  large  single  crystals  relatively 
easily. 

YIG  crystals  are  commercially  available  in  very  pure  form. 

1 I . B  Magnetic  Spin  Wave  Theory 

In  the  presence  of  an  applied  magnetic  field  (H),  a  magnetization 
field  (M)  is  induced  in  f errimagnetic  substances.  This  magnetization 
field  is  caused  by  the  alignment  of  individual  electron  spin  magnetic 
m.oments  and  is  effectively  a  volume  magnetic  dipole  moment.  If  a  static 
magnetic  field  is  applied,  then  both  magnetization  and  magnetic  fields 
are  present  inside  the  crystal.  The  magnetization  reaches  saturation 
when  essentially  all  the  electron  spin  moments  are  aligned. 

Magnetic  spin  waves  in  f errimagnetics  are  a  consequence  of  the 
existence  of  this  magnetization  field  and  its  incorporation  into  the 


constitutive  relations.  Maxwell's  Equations  for  a  sourceless,  materially 


isotropic  medium  (in  MKS  units)  are: 


VD  =  0 


VB  =  0 


V  X  H  =  BD/at 


V  X  E  =  -aB/at 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


The  constitutive  relations,  however,  are  defined  differently  from 
the  electromagnetic  free  space  case  in  that  (a)  the  magnetic  flux 
is  rrcportional  to  the  sum  of  the  magnetic  and  magnetization  fields 
(Eciiatior.  2.5),  and  (b.)  time  varying  magnetic  fields  interact  with  the 
magnetization  vector  as  torques,  inducing  a  time  dependence  in  the 
rr.,-;' r-e ;  i  zat  ion  (Equation  2.6  and  Figure  2.2): 


H  =  Pp(K  +  M) 


^  =  14. (M  *  Hgpp) 


=  H,,  -t-  h^,^  +  H, 


(2.5) 


(2.6) 


(2.7) 


~  2.8  X  10^  [Hz/Oe]:  gyroroagnetic  constant 


Rex  =  ^x’  R 

=  3  X  10" [m^]:  spin  exchange  constant 


^T***l**  f*rTPt 


M(c+dt ) 


M(t) 


Figure  (?.r) 


F'r'ocession  of  iragnetization  around  axis  of  static 
'  acnetic  field 


Equation  2.6  describes  the  interaction  between  the  internal  magnetic 


field  and  the  r.agnet  izat  ion .  The  presence  of  a  time  derivative  in  the 


constitutive  relation  causes  anisotropy  in  the  dispersion  relation  --  as 


will  be  seen  subsequently,  the  group  velocity  and  wave  vector  will  not  be 


collinear.  As  shown  by  Equation  2.8,  the  magnitude  of  M  is  constant  and 


the  cone  angle  (sec  Figure  2.2)  changes  as  the  driving  field  changes: 


8M  .  - 
~  •  M 
n  t 


Mr  =  0 


(2.8) 


Igpp,  the  effective  internal  magnetic  field  in  Equation  2.7,  can 


be  effectively  represented  as  containing  three  terms:  (a)  the  static 


applied  magnetic  field,  H^,  (b)  the  applied  RF  magnetic  field,  h^.^, 


and  (c)  the  exchange  interaction  field,  Hpy .  Additional  terms  are 


present  and  discussed  in  Appendix  A,  but  they  are  neglected  in  this 
approxiraat ion . 

If  it  is  assumed  that  (a)  the  DC  components  of  the  magnetization  and 
magnetic  fields  are  z  directed,  (b)  the  DC  components  are  much  larger 
than  the  time  varying  components,  and  (c)  the  time  varying  components  can 
be  modeled  as  plam 

M  =  MjjZ  +  fn 
H  =  HgZ  +  h 

(H„  is  the  z  directed  internal  magnetic  field),  and: 


e  waves,  then: 


M  =  M, 


(Mq  ■=  1780  [Oe]  in  YIG) 


(2.9) 


H  =  H, 


h(f ,t) 
m(r,t) 


h^eJ 

fBoe-J 


(a)t-k*r ) 
(a>t-k»r  ) 


(2.10) 


then  substitution  into  the  constitutive  relation  and  subsequent 
manipulation  (Equations  A.10-A.13)  yields  the  tensor  relations  between 
the  small  signal  magnetization  and  the  driving  magnetic  field.  The 
coefficient  matrix  is  the  Polder  Tensor: 


2  2 


2  2 


0(0^ 


"z"m 


2  2 


"z  = 


M  ^0  0 


o>  -  27tf 


(2.11) 


(2.12) 


The  Polder  Tensor  is  hermitian  and  the  medium  is  gyroraagnetic . 


Is?’ 


Applying  the  plane  wave  approximation  to  Equation  2.6  yields  the 


small  signal  constitutive  relation: 


jwm  =  yp  z  X  (o)  h  -  0)  m  - 


(2.13) 


Similarly,  applying  the  plane  wave  approximation  to  Maxwell's  Equations 
(A.10-A.19)  yields  relations  for  the  RF  magnetic  field  as  a  function  of 
the  small  signal  magnetization  field: 


-k(k«ra)  +  k^m 


(2.14) 


These  relations  define  the  characteristic  waves  of  the  medium  and  lead 
to  the  magnetostatic  approximation.  It  is  worth  noting  that,  for 
propagation  along  the  axis  of  magnetization,  the  h  and  in  fields  are 
transverse  and  circularly  polarized.  When  propagating  perpendicular  to 
the  axis  of  magnetization,  the  h  field  is  elliptically  polarized  with  a 
longitudinal  component  and  a  component  transverse  to  both  the  axis  of 
magnetization  and  the  direction  of  propagation.  From  Maxwell's  equations 
for  a  plane  wave  in  an  electrically  isotropic  medium  (A.15-A.16), 
however,  it  is  clear  that  the  e  and  b  fields  remain  transverse  to  the 
direction  of  propagation  [21]. 
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I I . C  The  Magnetostatic  Approximation 

Rather  than  substitute  Equation  2.14  into  Equation  2.13  to 
obtain  a  complete  dispersion  relation  (Equation  A. 21,  Figure  A.l),  a 
"magnetostatic"  approximation  can  be  applied  which  greatly  simplifies  the 
dispersion  relation  and  which  is  applicable  to  the  delay  line  problem. 

As  derived  in  Appendix  A  (Equation  A. 20),  the  small  signal  electric  field 
is  generally  proportional  to  1/k  and  a  magnetostatic  approximation  can  be 
made:  for  |k|  >>  k^  the  electric  field  approaches  zero,  the  curl  of  the 

time  varying  magnetic  field  approaches  zero,  h  can  be  approximated  as  the 
gradient  of  a  scalar  potential  field,  and  a  magnetostatic  approximation 
relating  the  time  varying  magnetic  and  magnetization  fields  is  available 
(Equations  A.22-A.28): 


r  -k(k*m) 
h  ~  - - 

k^ 


(2.15) 


Although  the  small  signal  electric  field  is  approximated  as  zero,  e  =  0, 
the  small  signal  energy  flow  is  still  accurately  described  by  a  small  but 
nonzero  Poynting  Vector,  e  x  h. 

Substituting  Equation  2.15  into  Equation  2.13  gives  the  dispersion 
relation  for  an  infinite  ferrite  under  the  magnetostatic  approximation: 


=  (0,^ 
s  i  n  (  0 )  = 


k  X  M 

|k||H| 


+  in  ( 0 ) ) 


(2.16) 


r-  I  I  ^ ^ j- 


This  dispersion  relation  (Figure  2.3)  is  central  to  the  modeling  of 


magnetoelastic  delay  lines. 

Incorporation  of  the  radial  boundary  conditions  for  a  YIG  rod 
produces  an  approximate  dispersion  relation  derived  by  Doane  [22] 
(Figure  2.4): 


(2.17) 


Jo(^R) 


k  +  k 

zero^^  order  Bessel  function 


This  equation  resembles  Equation  2.16  with  a  substitution  of  the 
arithmetic  mean  for  the  geometric  mean.  The  radial  component  of  k 
is  equal  to  the  ratio  of  the  zero*"^  Bessel  function  root  to  the  rod 
radius,  as  required  for  continuity  of  normal  b  and  tangential  h  at  the 
cylindrical  crystal  face.  Kittel's  original  derivation  of  the  dispersion 
relation  for  a  rod  demonstrated  the  radial  wavenumber  boundary  condition 
but  included  the  end  face  boundary  conditions  [16].  Doane's  derivation 
assumes  that  the  axial  wavelength  is  much  smaller  than  the  rod  length  and 
finds  the  same  radial  wavenumber  boundary  condition. 

It  is  of  interest  to  superimpose  the  dispersion  relation  which 
includes  the  radial  boundary  condition  on  the  magnetostatic  dispersion 
relation  for  an  infinite  medium  (Figure  2.5).  For  small  values  of  k  the 
wave  behaves  as  if  propagating  radially,  and  for  large  values  of  k  the 
waves  propagate  axially. 


Magnetostatic  dispersion  relation 


I - - - Kz  [Rod/M] 

Figure  (2.5)  -  Superimposed  magnetostatic  dispersion  relations:  (A)  Infinite  ferrite  - 
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II. D  Properties  of  Macnetostatic  Spin  Waves 


The  magnetoelastic  delay  line  is  made  possible  by  special  properties 
of  magnetostatic  and  magnetoelastic  waves,  including:  (a)  the  dependence 

of  both  phase  and  group  velocity  on  the  local  magnetic  field  strength  and 

« 

the  wavenumber,  and  (b)  the  existence  of  turning  and  crossover  points. 

The  following  sections  constitute  a  brief  summary  of  some  of  these 
properties  and  their  importance  in  delay  line  conception. 


II.D.l  Backward  Magnetostatic  and  Spin  Exchange  Wave  Behaviors 


The  dispersion  relations  of  Equations  2.16  and  2.17  both  allow  for 
dramatically  differing  wave  behaviors  as  a  function  of  the  wavenumber,  k. 
This  distinction  is  defined  for  "backward  magnetostatic  waves"  (referred 
to  in  the  literature  as  Backward  Magnetostatic  Volume  Waves  or  BMSVW) 
and  "spin  exchange  waves,"  although  both  are  manifestations  of  the 
magnetostatic  approximation  for  magnetic  spin  waves.  It  is  these 
strikingly  different  behaviors  that  allows  the  magnetoelastic  delay 
line  to  function. 

For  values  of  |k|  >>  k^  (to  maintain  the  magnetostatic 
approximation)  but  ^  (the  spin  exchange  contribution 

negligible),  the  dispersion  relation  in  Equation  2.17  can  be 
simplified  to: 


0)  +  - 
z  2 


k^  +  k^ 


(2.18) 


The  group  velocity  for  a  backward  magnetostatic  wave  propagating  in  this 
region  is  found  from  Hamilton's  equations  (see  Chapter  4): 


g  3k 


For  values  of  |k|  between  10^-10^  (rad/ra  *]  (Figure  2.6): 


^r  = 

|vg|  ~  10^-10^  [m/sec. 


(2.19) 


(2.20) 

(2.21) 


The  group  velocity  is  large  and  opposite  in  sign  to  the  axial  wavenumber. 
The  direction  of  propagation  of  energy  is  thus  opposite  to  the  direction 
of  phase  advance  and  the  wave  is  a  "backward  magnetostatic"  wave.  Also, 
the  components  of  group  velocity  are  not  directly  proportional  to  their 
corresponding  component  of  the  wave  vector.  The  group  velocity  and 
phase  vector  are  therefore  not  parallel  and  there  is  a  deflection  of 
propagation  direction  from  the  axis  of  phase  advance. 

For  values  of  |k|  >  5  x  10^  [rad/ra~'],  -  1  and  the  dipolar 

interaction  term  becomes  dominated  by  the  quantum  mechanical  exchange 
terra.  The  dispersion  relation  can  then  be  approximated: 


"  =  "z  ^  Vex^^ 


(2.22) 


This  behavior  is  delineated  in  Figure  2.6  as  the  "Spin  Exchange"  region. 

Applying  Hamilton's  equations  gives  the  components  of  group  velocity 
for  the  spin  exchange  dominated  wave: 


•  .•*  ."V  ..*•  -■.'h'"  .’"  l-N.  JV  -"o  Ji 


r<r*  •'w  Ti 


Magnetostatic  dispersion  relation  for  a  YIG  cyliner  (labeled) 
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(2.23) 

(2.24) 


The  group  velocity  components  are  directly  and  equally  proportional  to 
their  corresponding  components  of  the  vector.  The  group  velocity  is 
therefore  collinear  with  the  direction  of  phase  advance  and  the  wave 
is  a  "forward"  wave.  Even  for  very  large  values  of  the  wavenumber,  the 
spin  exchange  dominated  wave  has  low  group  velocity,  on  the  order  of 
10^-10^  [m/sec]. 


II. D. 2  Definition  of  Turning  Point 


The  transition  between  the  previously  described  "backward 
magnetostatic"  and  "spin  exchange"  wave  behaviors  is  critical  to 
understanding  the  operation  of  the  raagnetoelastic  delay  line.  Consider 
a  backward  magnetostatic  wave  propagating  axially  —  (a)  at  fixed 
frequency,  (b)  at  high  speed,  (c)  in  the  positive  direction,  and 
(d)  with  negative  wavenumber  --  toward  a  region  of  increasing  magnetic 
field  (Figure  2.7,  Point  A).  From  Equation  2.17  and  Figure  2.7,  the 
magnitude  of  the  (negative)  wavenumber  will  become  larger  and  the  wave 
will  move  "down"  the  dispersion  curve.  As  the  wave  continues  to 
propagate  towards  increasing  magnetic  field,  |k|  will  continue  to 
increase  until  the  wave  reaches  the  bottom  of  the  dispersion  curve 
(Figure  2.7,  Point  B).  At  this  point  the  wave  can  continue  to  propagate 
in  the  same  direction  only  with  imaginary  wavenumber.  This  point 
is  called  the  "Turning  Point"  and  it  is  here  that  the  "backward 
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magnetostatic"  behavior  of  the  wave  decreases  and  "spin  exchange" 
behavior  begins  to  dominate. 

The  increasing  spin  exchange  character  of  the  wave,  however, 
introduces  fundamental  changes  in  wave  behavior.  Since  the  spin  exchange 
wave  is  a  forward  wave,  the  group  velocity  changes  sign  and  becomes 

« 

negative  —  the  direction  of  energy  flow  actually  turns  around  and  the 
wave  propagates  back  toward  the  region  of  weaker  magnetic  field  (Figure 
2.7,  Point  C).  The  wavenumber,  however,  continues  to  be  negative  and 
increase  in  magnitude  —  the  "turning"  of  the  wave  is  unlike  a  reflection 
where  both  power  and  momentum  change  direction.  Also,  the  magnitude  of 
the  group  velocity  changes  dramatically,  from  backward  magnetostatic  wave 
values  of  10^  [m/sec]  to  spin  exchange  wave  magnitudes  of  lO’  [m/sec]. 

The  point  at  which  the  axial  group  velocity  equals  zero  as  it 
changes  from  positive  to  negative  is  the  "Turning  Point"  and  can  also  be 
defined  by  equal  and  opposite  spin  exchange  and  backward  magnetostatic 
contributions  to  group  velocity: 


2ii)w\  k 

n  ex  z 


tpl 


m„k^k^/k^ 

1/4 


=  0 


2\ 


ex 


(2.25) 

(2.26) 


Typically,  l^^pl  is  on  the  order  of  10^  [rad/m]  and,  from  Equations  2.12, 
2.17,  and  2.26,  the  magnetic  field  strength  is: 


(I) 


z 


(D 


(2.27) 
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ll. D. 3  Definition  of  Crossover  Point 

A  description  of  elastic  waves  and  their  coupling  to  magnetostatic 
spin  waves  is  contained  in  Appendix  A.  This  conversion  between  spin  and 
elastic  waves,  usually  with  high  efficiency,  occurs  at  the  "Crossover 
Point."  The  crossover  point  is  reached  when  the  wavenumber  reaches  the 
value : 


(2.28) 


v^^  =  3844  [m/sec]:  elastic  wave  velocity 


The  combined  magnetostatic  and  elastic  "magnetoelastic"  dispersion 
relation  is  shown  in  Figures  2.8  and  2.9. 


For  waves  at  microwave  frequencies  in  YIG,  values  of  wavenumber  at 
crossover  are  on  the  order  of  10^  [rad/m].  At  this  point,  the  character 
of  the  wave  changes  abruptly  from  spin  exchange  to  elastic  (or  vice 
versa).  As  described  in  Appendix  A. 3,  clastic  waves  are  forward  waves 
and  the  group  velocity  is  collinear  with  the  wave  vector. 


II. D. 4  Transition  Between  Turning  and  Crossover  Points 


From  the  previous  sections  it  is  possible  to  envision  a  magnetic 
spin  wave  whose  backward  magnetostatic  character  evolves  to  spin  exchange 
behavior  and  then  undergoes  transformation  to  an  elastic  wave.  These 
changes  would  occur  as  the  wave  propagates  into  a  turning  point,  reverses 
direction,  and  the  converts  into  an  elastic  wave  at  the  crossover  point. 
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The  cause  of  the  evolution  of  the  wave  would  be  the  presence  of  a 

spatially  inhomogenous  magnetic  field. 

The  striking  feaure  of  this  behavior  is  that  the  magnetic  spin  wave 

experiences  a  large  change  in  wavenumber  in  the  very  short  distance 

between  turning  point  and  crossover  point.  As  established  in  previous 

« 

sections,  at  the  turning  point: 


I  ktp I  =  10^  [rad/m] 


o/yyi^  [Oe] 


At  the  crossover  point,  |k  |  -  10^  [rad/m].  Applying  Equations  2.12, 

cp 


2.22,  and  2.28  gives  a  magnetic  field  strength  at  the  crossover  point: 


(2.29) 


For  a  nominal  frequency  of  4  GHz,  the  change  in  field  strength  between 
the  turning  and  crossover  points  is  a  very  small  22  Oe .  As  will  be  seen 
in  the  next  chapter,  derivatives  of  axial  field  profile  will  typically  be 
on  the  order  of  5  x  10^  [Oe/ra].  This  implies  that  the  magnetic  spin  wave 
undergoes  a  two  order  of  magnitude  increase  in  wavenumber  in  an  axial 


distance  on  the  order  of  0.1  [mm). 


Chapter  III 


THE  MAGNETOELASTIC  DELAY  LINE 


III. A  Introduction 


The  design  goal  for  the  magnetoelastic  delay  line  is  to  achieve  wave 
propagation  times  which  are  a  linear  function  of  frequency: 


TD(f)  =  D»f  +  T, 


(3.1) 


where  TD  is  the  time  delay,  D  is  the  dispersive  delay  coefficient,  and  T^^ 
is  some  fixed  delay.  This  linear  delay  function  has  numerous  practical 
appl icat ions . 

The  magnetoelastic  delay  works  by  converting  electromagnetic 
radiation  to  Long  wavelength  magnetostatic  waves  inside  a  YIG  crystal. 

The  YIG  crystal  has  a  nonuniform  internal  magnetic  field  designed  so  that 
the  magnetostatic  waves  will  evolve  into  spin  exchange  waves,  couple  with 
high  efficiency  to  short  wavelength  elastic  waves,  and  then  undergo  the 
exact  reverse  set  of  conversions,  ending  at  the  isolated  output  antenna. 
The  internal  magnetic  field  is  synthesized  to  produce  a  time  delay  which 
is  a  linear  function  of  frequency. 


III.B  One-Dimensional  Model  of  the  Maenetoelastic  Delay  Line 


The  operation  of  the  magnetoelastic  delay  line  can  be  understood  by 
following  the  path  and  evolution  of  a  wave  packet  of  energy  in  a  simple 
one-dimensional  model.  The  model  includes  the  following  features: 
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(a)  A  symmetric,  axially  directed,  concave  upward  magnetic  field  Hjj(z)Z 
is  present  (Figure  3.1a). 

(b)  A  constant,  axially  directed  magnetization  M^Z  is  present. 

(c)  The  medium  is  a  YIG  rod  with  the  dispersion  relation  of  Equation 
2.17  (Figure  2.4). 

(d)  Transducers  which  convert  electromagnetic  waves  to  spin  waves  are 
present  in  the  rod  but  do  not  interfere  with  the  propagation  of 
elastic  waves. 

The  delay  line  can  be  envisioned  as  operating  as  follows; 

-  Electromagnetic  energy  is  converted  to  a  long  wavelength,  high 
speed  backward  magnetostatic  wave  at  an  input  antenna  (Figure  3.1, 
Point  A).  Initially,  the  wave  has  <  0  and  propagates  in  the 
positive  z  direction  as  a  backward  wave  (Equation  2.20). 

-  As  the  wave  propagates  toward  greater  field  strength,  the  wavenumber 
increases  and  the  wave  slows  down  until  the  turning  point,  z^p 
(Figure  3.1,  Point  B),  is  reached.  At  z^.p,  the  wave  reverses 
direction  (Section  II. D)  and  becomes  a  forward  spin  exchange  wave 
traveling  in  the  -z  direction. 

-  The  wave  continues  to  travel  in  the  negative  z  direction  until,  as 
described  in  Section  II. D,  the  wavenumber  increases  in  magnitude 
sufficiently  to  permit  crossover  at  the  crossover  point,  z^p 
(Figure  3.1,  Point  C).  z^p  is  very  close  to  z^p  (within  0.1  [mm]). 

-  The  elastic  wave  continues  in  the  negative  direction,  travels 
through  the  center  of  the  rod,  and  then  the  entire  process  is 
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reversed  with  a  magnetostatic  wave  arriving  at  the  output  antenna 
(Figure  3.1,  Points  D-E) . 

The  evolution  of  the  magnetoelastic  wave  in  the  delay  line  is  shown  in 
Figure  3.2 


III.C  The  Delay  Function  and  Internal  Magnetic  Field 

Although  the  propagation  path  for  the  magnetoelastic  delay  line  has 
been  described,  the  means  to  produce  a  linearly  dispersive  time  delay  has 
not  been  defined.  This  section  surveys  the  rationale  for  synthesizing  a 
linear  delay  function  by  .'specifying  a  particular  on-axis  magnetic  field. 

The  total  time  for  a  wave  packet  to  travel  from  the  input  antenna 
to  the  output  antenna  in  the  manner  described  in  Section  III.B  is 
readily  estimated  and  is  a  very  good  approximation  to  the  actual 
three-dimensional  rod.  The  time  spent  in  propagation  can  be  broken  down 
into  three  portions:  (a)  backward  magnetostatic  wave j  (b)  spin  exchange 
wave;  and  (c)  elastic  wave.  The  backward  magnetostatic  wave,  however,  is 
very  fast  and  makes  a  negligible  contribution  to  the  overall  trip  time. 
Thus , 


TD(f ) 


T  +  T 

spin  exchange  wave  elastic  wave 


(3.2) 


Since  the  elastic  wave  is  insensitive  to  local  magnetic  field 
conditions  (except  at  crossover),  the  time  spent  as  an  elastic  wave  is 


simply: 
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The  time  spent  as  a  spin  exchange  wave  (whose  group  velocity  depends  upon 
the  local  magnetic  field  strength)  is  derived  in  Appendix  B.  The  result 


IS : 


ss 
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(3.4) 


As  shown  in  Chapter  2,  Equation  2.29: 
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and  the  total  time  spent  is: 


(3.5) 
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As  shown  in  Apptndix  B,  by  carefully  locating  the  crossover  point 

z  as  a  function  of  frequency,  the  linear  dispersion  can  be  obtained: 
c  p 


cp 
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(3.7) 


where  C  is  an  arbitrary,  nonpositive  constant, 


Generalizing  z^p  to  cover  all  values  within  the  rod  permits  complete 
derivation  of  the  required  one-dimensional  magnetic  field.  Synthesis  of 
this  on-axis  magnetic  field  produces  the  desired  linearly  dispersive 


function: 


(2|z|  -  Vg^T^)  +  '<2|z|  -  -  4v2^DC 


The  presence  of  an  arbitrary  constant  of  integration,  C,  is  a 
critical  factor.  Theoretically,  C  can  be  any  nonpositive  number  and  it 
determines  the  extent  of  nonlinearity  of  the  field  profile.  For  C  =  0, 
the  field  profile  is  a  linear  function  of  z.  As  C  becomes  m.ore  negative 
the  field  becomes  more  nonlinear.  Any  delay  function  therefore  has  a 
potentially  infinite  number  of  acceptable  field  profiles.  Optimizing  the 
selection  of  the  field  profile  constant  C  so  as  to  maximize  bandwidth  is 
one  objective  of  this  thesis. 


III.D  Insights  into  Delay  Line  Operation 


Not  all  the  magnetic  spin  wave  energy  arriving  at  the  output  antenna 
is  converted  back  into  electromagnetic  energy:  a  fraction  of  that  energy 
continues  to  propagate  until  it  reaches  the  input  antenna  and  repeats  the 
delay  process.  The  presence  of  second  (and  higher)  order  echoes  may 
limit  the  potential  bandwidth  of  the  delay  line.  For  example,  the 
presence  of  pulse  second  time  around  echoes  could  swamp  high  frequency 
leading  edge  transients  of  interest.  A  bandwidth  constraint  can  be 
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defined  by  requiring  that  the  second  order  echo  time  for  the  lowest 
frequency  of  interest  must  be  greater  than  the  delay  time  of  the  highest 
frequency  of  interest  (Figure  3.3); 


Figure  (3.3)  -  Echo  limited  bandwidth 


Dfhi  ^  <  2(Df^^  T^) 


(3.9) 


fy,,-  <  2f,  +  — 

hi  lo  D 


A  potential  maximum  bandwidth  can  be  readily  calculated. 

Also,  the  amount  of  time  spent  as  an  elastic  wave  and  as  a  spin  wave 
can  be  explicitly  solved: 


The  fixed  delay  term  T^^  is  wholly  incorporated  into  the  elastic  wave 
propagation  whereas  the  dispersive  propagation  time  Df  is  divided  equally 
between  the  elastic  and  spin  wave  regimes.  The  constant  of  integratiion 
term  C  adjustably  decreases  the  elastic  wave  propagation  time  while 
adding  equally  to  the  spin  wave  propagation  time. 

Since  the  spin  wave  propagation  time  is  explicitly  defined  in 
Equation  3.4  and  is  accounted  for  in  Equation  3.11,  a  simple  expression 
for  the  axial  derivative  of  the  magnetic  field  at  crossover  can  be 
obtained : 
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(3.12) 


The  magnetic  field  slope  at  crossover  is  independent  of  the  fixed  delay 
terra  T^^  and  is  inversely  proportional  to  C  and  D. 


1 1 1 . E  Focusing  Conditions 


Initial  work  in  raagnetostatic  spin  wave  delay  lines  found  that  the 
attenuation  of  the  signals  varied  strongly  with  the  magnetic  field 
profile's  characteristics  [23].  Theoretical  analysis  of  the  behavior  of 
spin  waves  between  z^p  and  z^p,  and  ray  tracing  simulations  demonstrated 
that  the  spin  waves  undergo  wavefront  focusing/defocusing  effects  which 
are  dependent  upon  the  local  magnetic  field.  Rays  which  defocus  show 
increased  attenuation  since  the  wavefront  disperses  through  the  rod. 
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Morgenthaler  has  analytically  solved  for  radial  position  and  the 


ratio  of  the  wavenumber  components  as  a  function  of  time  between  z^p  and 

z  [24].  His  analysis  shows  that  there  is  a  focusing  parameter  Q  which 
cp 

controls  whether  a  ray  will  (a)  focus,  (b)  defocus,  or  (c)  oscillate  in 

an  overfocused  condition.  Q  is  defined  by  the  magnetic  fields  at  the 

« 

crossover  point,  z_„: 
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(3.13) 


For  Q  <  0,  the  wave  is  defocusing;  Q  >  1  is  overdamped;  and  1  >  Q  >  0  is 
a  focusing  condition. 

The  one  dimensional  delay  line  approximations  in  Equation  3.12  can 
be  applied  to  the  definition  of  Q.  The  result  is  an  explicit  expression 
for  Q  as  a  function  of  the  selected  linear  dispersion  function,  the 
frequency,  and  the  field  profile  constant: 


Q  = 


3a 

(a  +  1)2 


4a  r 
a  +  1  2C 
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(3.14) 


For  the  linear  delay  line,  Q  can  range  between  1  and  An  overfocusing 

condition  is  not  possible  with  these  types  of  profiles. 

Since  Q  is  not  defined  for  negative  magnetic  field  gradients,  it  is 
not  generally  applicable  to  the  complete  delay  line  problem.  It  does. 
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however,  provide  a  design  criterion  which  correlates  with  bandwidth  (as 
will  be  shown).  It  is  worth  noting,  however,  that  D  must  be  greater  than 
zero  for  a  focusing  profile  (Q  >  0)  to  be  achieved. 


III.F  Synthesis  of  the  Magnetic  Field 


Section  III.C  derives  the  general  formula  for  the  on-axis  ma^uetic 
field  for  any  given  linear  delay  function.  Appendix  C  reviews 
Morgenthaler ' s  approach  to  synthesis  of  any  on-axis  field  in  a  rod  or 
disk.  The  field  is  created  by  manufacturing  pole  pieces  whose  surfaces 
are  continuous  with  the  equipotential  which  would  be  generated  by  the 


specified  field  (Figures  3.4,  3.5,  and  3.6).  The  pieces  are  spherically 
symmetric,  made  from  ferrous  material,  and  energized  in  an  external 
magnetic  field. 


Figure  (3.5)  -  Pole  piece  design 


Computer  generated  equipoten 
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Chapter  IV 
RAY  TRACING 


IV. A  Introduction 


Ray  tracing  is  the  primary  technique  used  herein  for  analysis  of 
the  isolated  two  port  magnetoelastic  delay  line.  Ray  tracing  has, 
of  course,  been  the  primary  tool  in  optical  analysis  for  over  a  century 
and  predates  Maxwell's  equations.  Ray  tracing  in  anisotropic  media, 
however,  is  somewhat  more  complex  and  relies  heavily  upon  some  additional 
assumptions.  This  chapter  will  explore  the  methodology  of  ray  tracing, 
the  underlying  assumptions,  and  the  computer  program  developed  for  this 
thesis . 


IV. B  Theory  of  Ray  Tracing 

Ray  tracing  is  based  on  the  idea  that  if  an  electromagnetic  wave 
propagating  through  a  medium  meets  certain  conditions,  it  can  be  treated 
as  a  particle  which  obeys  Newtonian  laws  of  motion,  i.e.,  the  particle 
motion  can  be  approximated  by  a  fixed  velocity  over  small  time  intervals. 
This  approximation  allows  the  tracing  of  ray  paths  which  are  the 
"geodesics  for  the  metric  where  the  optical  length  ...  is  a  measure  of 
the  arc  ..."  [25].  The  local  nature  of  the  geodesic  depends  upon  the 
refractive  index  of  the  material. 

Geometrical  optics  has  been  shown  to  be  a  valid  way  of  describing 
the  propagation  of  magnetoelastic  waves  in  nonuniform  magnetic  fields 
[20,26].  Morgenthaler  has  demonstrated  that  a  particle  approximation  to 
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magnetoelastic  wave  propagation  in  magnetically  saturated  ferro-  and 
ferrimagnets  is  a  natural  consequence  of  small  signal  energy  and  momentum 
conservation  laws  [20].  His  proof,  which  incorporates  the  terms  in  the 
complete  magnetoelastic  wave  description  given  in  Equations  A.5-A.7, 
requires  only  that  wavenumber  be  large  compared  with  both  the  normalized 
spatial  field  derivatives  and  the  rod  dimensions: 


H(r,2) 


(4.1) 


|k  I  >>  radius,  length 


(4.2) 


Given  these  conditions,  the  spatial  dependence  of  the  small  signal  fields 
can  be  represented  by: 


:r,t)  =  s  ej(^t-Ik»dr) 


(4.3) 


and  a  "quasiparticle"  description,  based  upon  a  linear  sum  of  plane 
waves,  is  acceptable. 

Morgenthaler  applied  the  above  assumptions  to  the  magnetoelastic 
wave  small  signal  power-energy  and  stress-momentum  conservation  theorems 
and,  after  taking  a  time  average  over  the  period,  obtained  the  Hamilton 
Equations  of  optics: 


dr  -  I 
dt  "  ^k‘^'r,t 


(4.4) 
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Equations  4.4  and  4.5  form  the  basis  for  ray  tracing  techniques. 
Equation  4.4  is  the  classical  group  velocity  relation  of  optics. 
Equation  4.5  gives  the  time  evolution  of  the  wavenumber  as  a  function 
of  the  spatial  derivative  of  the  frequency.  Ray  paths  are  found  by 
integrating  Equations  4.4  and  4.5;  although  most  ray  paths  are  actually 
found  from  numerical  integration  over  successive  small  time  increments, 
closed  form  solutions  are  sometimes  available. 

The  time  derivative  of  the  wavenumber  (the  group  velocity  solution 
has  already  been  presented  in  Section  II. D)  can  be  found  from  Equations 
2.18  and  2.22: 


dt 


z  + 


aH  _  n 
3r  ^  ^ 


(4.6) 


The  solutions  for  both  backward  magnetostatic  and  spin  exchange  wave 
behaviors  are  approximately  the  same. 

The  magnetic  fields  inside  the  rod  are  found  from  Taylor  series 
expansions  of  the  on-axis  field  (Appendix  C  and  Figures  4.1  and  4.2). 
Kear  the  central  axis,  only  the  highest  order  component  of  each  spatial 
derivative  needs  to  be  considered: 


|H(r,z) I 
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From  Equation  3.12,  the  axial  derivative  of  the  on-axis  magnetic  field  at 
the  crossover  point  can  be  approximated  and  substituted  into  both 
equations ; 
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The  spatial  derivatives  throughout  the  YIG  rod  remain  constant  to  within 
an  order  of  magnitude,  the  rate  of  time  evolution  of  the  axial  wavenumber 
is  greater  than  that  of  the  radial  wavenumber,  i.e.,  |dk^/dt|  >> 

|dk^/dt|,  and  the  magnitude  of  k^  becomes  much  larger  than  the  magnitude 
of  k^.  It  is  in  this  manner  that  spatial  inhomogeneity  of  the  magnetic 
field  controls  the  time  evolution  of  wavenumber. 


IV . C  Dynamic  Range  and  Saturation 

Ray  tracing  can  also  be  used  to  estimate  power  handling  capability 
up  to  levels  where  saturation  effects  produce  nonlinearities.  Since  a 
finite  area  contains  a  limited  number  of  processing  atoms,  the  magnetic 
spin  wave  energy  density  cannot  exceed  the  maximum  level  called 
saturation.  Above  this  level,  higher  order  processes  are  present  which 
produce  nonlinear  effects  [27). 

For  any  plane  wave  based  on  Maxwell's  Equations,  the  power  flow 


through  a  unit  normal  cross-sectional  area  (Figure  4,3)  is  proportional 
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(4.10) 


Figure  (4.3)  -  Power  flow  from  ray  tracing 


=  intensity 
=  cross-sectional  area 
=  power 

=  Idz  cr  VglAl^dz 
=  field  magnitude 


Ray  tracing  Jends  it^plf  well  to  estimating  relative  power  saturation 
points  through  application  of  Equation  4.10,  since  the  distance  between 
rays  can  be  related  to  cross-sectional  area.  Magnetoelastic  waves  are 
relatively  low  loss  and  the  field  amplitude  can  be  assumed  to  be 
constant.  The  power  can  therefore  be  considered  proportional  to  the 
product  of  cross-sectional  area  and  group  velocity.  Saturation  occurs 
in  the  ray  paths  where  Equation  4.10  is  a  minimum. 


IV. D  Considerations  in  Ray  Tracinc 


Two  important  factors  in  ray  tracing  are:  (a)  the  causes  of 
ray  path  bending,  and  (b)  the  behavior  of  rays  in  isotropic  versus 
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anisotropic  media.  Brief  consideration  of  these  issues  is  worthwhile  for 
the  insight  into  ray  tracing  in  a  gyrotropic  medium. 

There  are  two  primary  causes  for  ray  paths  departing  from  straight 
line  propagation  in  gyromagnetic  media.  First,  the  dispersion  relation 
is  anisotropic  --  the  group  velocity  and  wave  vectors  are  not  collinear. 
This  deflection  phenomenon  is  illustrated  by  considering  the  wave  surface 
for  a  gyrotropic  medium.  The  wave  surface  in  Figure  4.4  is  defined  for 
the  magnetostatic  dispersion  relation  of  Equation  2.17  at  fixed  frequency 
and  magnetic  field  strength.  It  can  be  proven  that  the  group  velocity  is 
always  perpendicular  to  the  wave  surface  [28]  but  is  not  necessarily 
parallel  to  the  wave  vector: 


Kz 


Figure  (4.4)  -  Magnetic  Wave  Surface  and  Anisotropy 

Clearly,  a  small  shift  in  the  displayed  angle  q  of  the  wave  vector 
produces  a  dramatic  deflection  in  the  direction  of  the  group  velocity  and 
an  effective  bending  of  the  ray. 

The  magnetic  inhomogeneity  of  the  medium  also  produces  ray 
refraction.  A  ray  propagating  along  the  central  axis  of  the  rod 
experiences  a  greater  magnetic  field  than  a  ray  propagating  off-axis. 

For  fixed  frequency,  rays  propagating  between  the  turning  point  and 


crossover  point  may  have  dramatically  different  phase  velocities  but 
similar  group  velocities,  as  shown  in  Figure  4.5.  The  result  is  that  the 
phase  front  tends  to  refract  and  ray  paths  bend. 

The  use  of  wave  surfaces  as  shown  in  Figure  4.4  is  an  extremely 

useful  means  for  visualizing  how  ray  paths  evolve.  Before  going  on  to 

« 

analyze  the  magnetoelastic  delay  line,  however,  it  is  worthwhile  to 
examine  a  few  simple  cases. 


Figure  (4.6)  -  ''/ave  Surface  for  Electromagnetic  Wave  in 
isotropic  Medium 

The  case  shown  in  Figure  4.6  is  the  classical  optics  case, 
propagation  of  light  in  a  uniform  medium.  Because  the  medium  is 
homogenous  and  the  dispersion  relation  is  isotropic  (the  constitutive 
relations  contain  no  time  derivatives),  the  group  and  phase  velocities 
are  collinear. 
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FigurG  (4.7)  - 


Electromagnetic  Radiation  in  a  Waveguide  with  Spatially 
Dependent  Permittivity 


The  case  shown  in  Figure  4,7  illustrates  the  behavior  of 
electromagnetic  radiation  in  a  waveguide  filled  with  a  substance 
wiiose  permittivity  increases  with  x.  The  transverse  component  of 
the  wavenumber  is  fixed  by  the  boundary  condition.  The  longitudinal 
component  decreases  until  it  reaches  a  zero  cutoff;  at  that  point  the 
wave  is  reflected  although  some  energy  continues  to  propagate  as  a 
decaying  exponential.  Spatial  dependence  of  permittivity  is  analogous 
to  the  spatial  dependence  of  susceptibility  in  the  magnetoelastic  delay 
line.  The  wave  surface  can  be  imagined  as  a  circle  whose  radius  shrinks 
with  increasing  permittivity. 


source 


0 


iou  0  '4.-)  -  The  Luneberg  Lens 


Figure  4.8  shows  the  path  of  rays  in  a  Luneberg  Lens,  commonly  used 
in  antenna  systems  where  aperture  blockage  is  unacceptable  [29].  The 
Luneberg  Lens  focuses  a  plane  wave  onto  a  surface  point  through  the  use 
cf  a  spatially  varying  permittivity: 


.  (r)  =  2  -  (  ^  )^ 


(4.11 
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The  wave  surface  for  magnetoelastic  waves  (fixed  frequency,  field 
strength)  is  shown  in  Figure  4.9.  The  radial  portions  of  the  surface 
correspond  to  magnetostatic  behavior,  the  circular  end  portions  to  spin 
exchange  wave  behavior.  The  area  contained  inside  the  surface  decreases 
as  -*  CD.  Clearly,  backward  magnetostatic  wave  group  velocities 
experience  strong  deflection  from  their  wave  vectors;  spin  exchange 
waves  have  collinear  group  velocity  and  wave  vectors. 


The  wave  surface  does  not  give  any  indication  of  group  velocity 
magnitude.  This,  however,  is  easily  visualized  from  the  m-k  dispersion 
diagram,  since  the  group  velocity  is  the  slope  of  the  dispersion  curve. 

Table  4.1  lists  several  key  quantitites  needed  to  visualize  the  ray 


tracing  of  the  magnetoelastic  delay  line: 
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Table  4.1  Summary  of  Backward  Magnetostatic  Wave,  Spin  Exchange  Wave, 
and  Elastic  Wave  Relations 

Initial  value  of  k^.  is  defined  by  Kittel's  analysis  of  magnetostatic 
propagation  in  rods  [16],  An  initial  value  of  k^  is  found  by  solving  the 
magnetostatic  (without  exchange)  dispersion  relation  in  Equation  2.18. 
Froms  Equations  4.8  and  4.9,  it  is  apparent  that  k^  will  evolve  much  more 
quickly  than  k^,. 

Figures  4.10a-d  show,  in  parallel,  the  evolution  of  the  ray  path 
and  the  corresponding  visualization  from  the  wave  surface.  Since  dkj./dt 
is  always  positive  and  k^.  is  initially  negative,  the  evolution  is  not 
symmetric.  The  evolution  of  k  ,  however,  is  roughly  symmetric  as  dk  /dt 


changes  sign  with  the  sign  of  z. 


A.  Startin 


backward  ma 


The  ray  tracing  computer  program  iteratively  tracks  out  the  ray  path 
using  the  techniques  described  above.  The  dispersion  relation  used  is 
that  of  Equation  2.16,  r.agnctostat ic  waves  in  an  infinite  ferrite.  Use 
of  this  dispersion  relation  is  acceptable  because  the  initial  conditions 
guarantee  that  the  wavelength  is  virtually  always  much  smaller  than  the 
rod  dimensions  and  the  medium  appears  infinite  to  the  wave.  Borgeaud's 
analysis  [3]  confirms  the  acceptability  of  (a)  applying  ray  tracing 
technqiues  to  relatively  long  wavelength  raagnetoelastic  waves  in  the 
early  stages  and  (b)  using  the  first  root  of  the  zero^^  order  Bessel 
function  as  the  initial  value  of  k^.. 

The  Hamilton  equations  are  applied  to  this  dispersion  relation 
without  any  simplifying  assumptions.  The  conditions  for  application  are 
that  (a)  frequency  is  constant,  and  (b)  only  H  and  M  are  functions  of 
position.  In  the  latter  case,  fl  and  R  are  defined  anywhere  in  the  rod  by 
the  method  described  in  Appendix  C.  The  ray  tracing  is  implemented  as 
shown  in  Figure  A.  11.  At  Point  1,  the  position,  wave  vector,  and  group 
velocity  are  known.  For  a  time  increment  on  the  order  of  nanoseconds,  Ar 
and  Ak  are  calculated  from  Hamilton's  equations.  The  new  position,  wave 
vector,  and  group  velocity  are  then  found. 


TD(f)  -  400f  -  100  [noGc]  Q  =  0.27 


Figure  (4.12)  -  Computer  generated  ray  tracing  for  bandwidth  analysis 
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Chapter  V 


RESULTS  AND  CONCLUSIONS 


V.A  Introduction 


Ray  tracing  analysis  provides  several  insights  into  the  operation 
and  optimization  of  the  two-port  "buried-antenna"  magnetoelastic  delay 
line.  As  a  result,  bandwidth  and  dynamic  range  can  be  optimized,  and  the 
effects  of  several  different  material  parameters  are  explained.  A  simple 
predictive  model  is  developed. 


V'.B  Summary  of  Previous  Work 


Borgeaud  [3]  built  and  tested  a  two-port  "buried-antenna" 
.-■nagnetoelastic  delay  line  (Figure  3.4)  designed  to  the  following 
parameters : 


TD(f)  =  500  f  -  200  [nsec. 


C  =  -1000 


Rod  radius  =  1.5  [mm] 


Rod  length  =  5  [mm] 

Axial  antenna  locations  =  +/-  1.0  [mm] 


The  value  of  C  =  -1000  was  chosen  to  demonstrate  that  a  nonlinear  field 
profile  would  produce  a  linearly  dispersive  delay  function.  His  results 
can  be  summarized  as  follows: 


SSI 


Range  of  detectable  signal  “  3. 1-3.6  [GHz] 

Detectable  signal  bandwidth  «  500  [MHz] 

Range  of  linear  dispersion  =  3. 2-3. 5  [MHz] 

Linearly  dispersive  signal  bandwidth  =  300  [MHz] 

Actual  delay  function  implemented  =  500  f  -  500  [nsec] 


Borgeaud  found  that  the  range  of  frequencies  which  propagate  is 
determined  by  the  antenna  location.  Only  those  frequencies  falling 
within  the  range  defined  below  can  propagate: 


(5.1) 


This  region  is  the  spin  wave  manifold  and  is  bounded  by  the  dispersion 
relation  curves  shown  in  Figure  5.1.  In  this  case  is  defined  by  the 
local  field,  at  the  antenna.  Borgeaud's  selection  of  antennas  at 

+/-  1.0  [ram]  explicitly  defined  the  available  range  of  frequencies  for 
the  given  magnetic  field  profile.  For  &>  -*  the  frequency  of 

propagation  is  already  at  the  turning  point  and  the  rays  propagate 
immediately  toward  regions  of  lower  magnetic  field  (Figure  5.2a);  as 
frequencies  increase  the  rays  spend  more  time  as  a  backward  magnetostatic 
wave  and  propagate  closer  to  the  z-axis  (Figure  5.2b).  Rays  which 
propagate  away  from  the  z-axis,  however,  show  large  deviations  from  the 
desired  linear  time  response  (Figure  5.2a,b). 


•  -* 
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Magnetostat 


V. C  Bandwidth  Analysis 

In  order  to  analyze  the  bandwidth  obtained  by  Borgeaud,  it 
was  necessary  to  create  a  computer  program  which  would  conduct  the 
appropriate  ray  tracing.  This  program,  described  in  Section  IV. F,  meets 
design  goals  of  (a)  a  lOX  increase  in  computation  speed  over  earlier  ray 
tracing  routines,  and  (b)  automatic  bandwidth  analysis.  Overall  program 
speed  is  improved  by  hosting  on  a  faster  microcomputer,  writing 
computationally  intensive  routines  in  C  language  (instead  of  Basic), 
and  optimizing  code  structure. 

For  the  range  explored  by  Borgeaud  (a  bottom  frequency  of 
3.1  [GHz]),  it  is  apparent  that  larger  bandwidths  are  obtainable  by 
changing  the  value  of  the  field  profile  constant,  C.  Figure  5.3  shows 
how  bandwidth  varies  as  a  function  of  C  for  several  different  delay 
functions  and  the  same  bottom  frequency.  The  antenna  location  explicitly 
defines  the  resonant  frequency  for  the  given  frequency  range.  A  relation 
riving  bandwidth  as  a  function  of  C  and  D,  however,  is  not  readily 
apparent . 

Although  the  choice  of  D,  the  dispersive  delay,  has  a  marked  effect 
on  bandwidth,  the  choice  of  the  fixed  delay  Tp  has  no  effect  on  bandwidth 
'Figures  f.4a  and  5.4bl.  This  implies  that  the  bandwidth  is  independent 
of  the  location  of  the  crossover  point  (Equation  3.7)  —  an  intuitively 
appealing  result  since  the  elastic  wave  portion  of  a  ray  path  tends  to  be 
roughly  parallel  to  the  central  axis.  The  elastic  wave  portion  of  the 
propagation  does  not  constrain  bandwidth  and  the  fixed  delay  T^^  is 


contained  entirely  in  the  elastic  wave  propagation. 


SAMPLE  DATA 


^  s 


Figure  5.3  illustrates  that,  for  each  value  of  dispersive  delay  D, 


the  bandwidth  reaches  a  maximum  for  some  value  (Figure  5.5a-c).  For 

values  of  C  more  negative  than  higher  frequency  rays  collide  with 

the  z-axis  and  never  reach  the  output  antenna  (Figure  5.6).  For  values 
of  C  greater  than  C^^x’  between  the  se.:.ond  turning  point  and  the 

output  antenna  experience  abrupt,  unusual  changes  in  direction  (Figure 
5.7).  The  former  case  is  clear  from  ray  tracing;  the  latter  effect 
appears  due  to  the  interaction  between  k  (which  evolves  asymmetrically) 
and  the  radial  components  of  M. 

It  was  suggested  by  Morgenthaler  [30]  that  bandwidth  maxima  might 
be  dependent  upon  Q,  the  focusing  parameter  discussed  in  Section  III.E. 
Figure  5.8  shows  this  to  be  the  case,  with  maximum  values  of  bandwidth 
consistently  occurring  for  0.4  <  Q  <  0.5.  Rays  with  a  value  of  Q  greater 
than  0.5  tend  to  show  the  abnormal  changes  in  direction  between  the 
second  turning  point  and  the  output  antenna.  Frequencies  whose  rays  have 
values  of  Q  less  than  0.5  obey  a  regular  trend. 


V . D  A  Simple  Predictive  Model 


The  independence  of  bandwidth  from  T^  leads  to  the  following 
!.y;  ' .  thes  is :  Bandwidth  is  a  function  of  dH/dz.  Tiie  inference  is  made 
ije^'acse  tiie  slope  of  the  magnetic  field  at  crossover,  dh(r  =  0,z)/dz  at 
z  =  Z|,p,  is  also  independent  of  T„  (Equation  3.12): 
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Figure  (5.6)  -  Ray  propagates  directly  from  antenna  to 
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The  derivative  of  magnetic  field  in  Equation  5.2  is  a  function  of  D  and 
C,  and  it  remains  constant  to  within  an  order  of  magnitude  throughout  the 
rod.  A  Log-Log  plot  of  bandwidth  versus  dH/dz  (excluding  points  with 
Q  >  0.5)  shows  a  straight  line  relationship  (Figure  5.9a)  with  a  slope  of 
2.3  (from  least  mean  square  line  fitting).  A  linear  plot  of  bandwidth 

t 

versus  (dH/dz)^'^  gives  a  percentage  standard  deviation  of  10%,  small 
enough  for  predictive  use  (Figure  5.9b). 

The  dependence  of  bandwidth  on  the  slope  of  the  magnetic  field 
profile  seems  reasonable  because  dk^/dt  and  dk^/dt  are  proportional  to, 
respectively,  dH/dz  and  dH/dr  (Equation  4.6): 

(5.3) 

(5.4) 
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Recall  also  that  dH/dr,  from  Equation  4.7,  can  be  expressed  as  a  function 
of  the  axial  field  derivative. 
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As  mentioned  above,  the  higher  the  frequencies  radiated  from  a  fixed 
antenna  location,  the  closer  the  rays  propagate  to  the  z  axis.  Above  the 
upper  edge  of  the  frequency  band,  rays  collide  with  the  z  axis.  The 
controlling  factor  for  the  backward  magnetostatic  wave  is  thus  the 
evolution  of  v^  (Equation  2.21): 


A.  LOG-LOG  plot  of  bandwidth  vs  dH/dz 


Figure  (5.9)  -  Bandwidth  as  a  function  of  dH/dz 
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(5.6) 


All  other  things  being  equal,  the  faster  the  magnitude  of  approaches 
zero  and/or  becomes  positive,  the  less  likely  the  ray  is  to  collide  with 
the  z-axis.  For  z  >  0:  and  k^  are  initially  negative;  dk^/dt  <  0; 

dk^/dt  >  0;  and  iarger  magnitudes  of  dkj./dt  and  dk^/dt  improve  the 
likelihood  of  a  ray  reaching  the  turning  point  and,  ultimately,  reaching 
the  output  antenna. 

Borgeaud's  assumption  that  rays  reaching  the  z-axis  propagate  no 
further  is  an  intuitive  consideration  of  energy  flow.  Rays  reaching  the 
z-axis  from  opposite  ends  of  the  same  antenna  tend  to  intersect  there 
and  ray  paths  can  no  longer  be  defined.  Ray  tracing  analysis  of  rays 
reaching  the  z-axis  confirms  the  assumption  on  different  grounds.  Those 
rays  which  have  v^  <  0  at  the  z  axis  do  not  appear  to  ever  leave  the 
axis.  Figure  5.10  shows  the  behavior  of  such  a  ray  --  v^  remains 
negative  and  the  ray  oscillates  about  the  z  axis  with  increased 
attenuation . 

Experimenting  with  bandwidth  data  and  Equations  5.2  and  5.5  leads  to 
the  following  proportionality: 
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This  relation  preserves  the  order  of  the  proportionality  between 
bandwidth  and  dH/dz,  and  it  allows  prediction  of  bandwidth  as  a  function 
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gure  (5.10)  -  Detail  of  ray  reaching  z-axis  and  osci 


of  several  other  variables:  (a)  Mg,  the  saturation  magnetization; 

(b)  the  elastic  wave  velocity;  and  (c)  yjig,  the  gyromagnetic 

constant.  Figures  5.11a-c  show  ray  tracing  bandwidth  data  which  are  in 
accordance  with  expected  values. 

As  a  result  of  this  work,  optimizing  bandwidth  for  a  given  delay 
function  becomes  a  trivial  calculation.  Recasting  Equation  3.14  to  give 
C,  the  field  profile  constant,  as  a  function  of  Q  allows  immediate 
calculation  for  Q  =  0.5; 


C 
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.  Q(a  +  1)^  -  3a  -  4a(a  +  1)  . 


(5.8) 


This  value  optimizes  bandwidth. 

V . E  Dynamic  Range 

A  second  goal  of  this  thesis  is  to  determine  if  the  dynamic  range 
of  the  magnetoelastic  delay  line  can  be  optimized  as  well.  Dynamic  range 
is  defined  here  in  the  sense  of  a  saturation  point  since  no  study  has 
yet  been  made  of  noise  floor  and  Minimal  Discernible  Signal  (MDS).  As 
discussed  in  Section  IV. C,  power  flow  can  be  viewed  as  a  function  of 
the  distance  between  rays  emanating  from  opposite  ends  of  the  transmit 
antenna.  The  power  handling  capability  is  proportional  to  the  product  of 
cross-sectional  area  between  rays  and  the  group  velociy,  assuming  the 
fields  are  of  constant  amplitude. 

Figure  5.12  shows  the  paths  of  rays  for  the  same  frequency  and  delay 
functions  but  in  magnetic  fields  with  different  values  of  C.  The 


distance  between  the  rays  decreases  as  C  approaches  zero  and  the 


distances  remain  relatively  constant  past  the  turning  point.  Since  the 
group  velocity  at  the  turning  points  is  very  small  and  similar  for  both 
sets  of  ray  paths  (~  1  [m/sec]),  dynamic  range  should  be  improved  for  the 
magnitude  of  C  increasing.  This  is  opposite  to  the  results  for 
bandwidth . 


V . F  Consequences  of  Field  Asymmetry 

The  development  of  a  comprehensive  ray  tracing  program  has  allowed 
several  other  aspects  of  magnetoelastic  delay  line  behavior  to  be 
explored.  One  issue  of  interest  is  the  asymmetry  between  ray  paths 
in  the  delay  line.  In  defining  the  delay  line,  the  magnetic  and 
magnetization  fields  are  defined  as  vectors  whose  axial  components  are 
directed  in  the  +z  direction.  The  results  of  Sections  V.B-C,  however, 
are  only  for  rays  injected  at  the  antenna  in  +z  half  plane  and  whose 
axial  component  of  the  elastic  wave  velocity  vector  is  negatively 
directed . 

The  absence  of  axial  magnetic  and  magnetization  field  symmetry  about 
z  =  0  implies  that  ray  paths  will  not  be  symmetric.  Likewise,  the  Taylor 
series  expansion  of  the  on-axis  magnetic  and  magnetization  fields  through 
the  rod  gives  radial  field  components  which  are  odd  with  respect  to 
z  =  0.  These  conditions  cause  differences  in  (a)  the  initial  conditions 
of  k^  and  (b)  the  values  of  k  x  B  and  k  •  fi  in  the  different  half-pJ.anes . 

Comparing  bandwidth  ray  tracing  for  rays  starting  at  the  +z  and  -z 
antennas  under  otherwise  identical  conditions  shows  that  the  behaviors 
are  very  similar,  if  not  identical.  Intuitively  again,  this  is  appealing 
since  dk^/dt  and  dkj./dt  are  roughly  symmetric  about  z  =  0  and  the  change 


in  Initial  value  of  is  small  compared  to  the  magnitude  of  the  total 
change.  The  changes  in  the  vector  products  have  minimal  impact  on 
both  the  group  velocity  and  dk/dt.  This  result  has  strong  practical 
implications:  the  delay  line  can  be  used  in  either  a  forward  or  backward 

sense  with  respect  to  the  direction  of  magnetic  field.  Choice  of  antenna 
terminal  has  no  apparent  effect. 

V . G  Comparison  to  Single-Ended  Delay  Line 

Strauss  type  delay  lines  with  a  synthesized  magnetic  field  can  also 
be  analyzed  by  the  ray  tracing  program.  In  fact,  the  differences  in 
computer  analysis  of  the  two-port  "buried-antenna"  and  two-port,  single- 
ended  Strauss  type  lines  are  trivial;  the  latter  merely  fixes  the  antenna 
at  2  =  0  and  requires  elastic  wave  reflection  at  z  =  0 .  The  delay  line 
built  by  Platzker  and  Morgenthaler  [2]  can  be  considered  in  this  manner 
and  good  correlation  with  experimental  results  is  found. 

Ray  tracing  demonstrates  that  the  two-port  "buried-antenna" 
design  provides  an  additional  degree  of  freedom  in  design  as  well  as 
substantially  larger  bandwidths.  The  degree  of  freedom  results  from 
being  able  to  locate  the  antenna  so  that  the  local  magnetic  field  is  at 
the  resonant  frequency.  In  the  Strauss  type  line,  the  antenna  is  fixed 
at  2  =  0,  and  T^  and  C  cannot  be  set  independently. 

The  lower  bandwidths  evidenced  for  the  Strauss  type  delay  line  are  a 
resvilt  of  waves  at  frequencies  near  the  bottom  of  the  spin  wave  manifold 
not  reaching  a  crossover  point  (Figure  5.13a).  Absent  this  low  frequency 


range  effect,  the  bandwidths  of  the  "buried-antenna"  and  Strauss  type 


delay  lines  would  be  the  same.  As  described  in  Section  V.A,  spin  waves 


( 
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at  or  near  the  local  resonant  frequency  in  the  two-port  "buried-antenna" 
line  are  immediately  at  the  turning  point  and  propagate  towards  regions 
of  lower  magnetic  field  (Figure  5.13b).  In  the  Strauss  line,  however, 
there  is  no  region  of  lower  magnetic  field  since  the  antenna  is  at  z  =  0. 
Waves  at  or  near  this  frequency  cannot  reach  a  crossover  point  and  simply 
return, to  the  output  antenna  as  a  magnetostatic  wave. 

Analysis  of  the  Strauss  type  delay  line  also  suggests  that  rays 
which  originate  near  the  center  of  the  endface  are  less  likely  to 
propagate  than  rays  starting  at  greater  radial  distance.  As  with  the 
two-port  "bur i ed-anter.na"  delay  line,  the  radial  velocity  of  a  ray  must 
reach  zero  before  the  ray  collides  with  the  z  axis.  Rays  starting  at 
greater  radial  distance  have  greater  likeihood  of  reaching  a  turning 
point  (Figures  5.14a  and  5.14b).  Rays  returning  from  reflection, 
however,  reach  the  endface  at  lesser  radial  distances  (Figure  5.14b). 

This  suggests  that  endface  antenna  designs  siiould  be  optimized  for 
maximum  coupling  at  different  radial  distances. 
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Chapter  VI 

CONCLUSIONS  AND  RECOMMENDED  FUTURE  WORK 

A  ray  tracing  analysis  of  "buried-antenna"  and,  to  a  lesser  extent, 
Strauss  type  linearly  dispersive,  two-port  magnetoelastic  delay  lines  has 
been  performed.  It  has  been  found  that  either  operating  bandwidth  or 
dynamic  range  can  be  maximized  with  respect  to  the  choice  of  internal 
magnetic  field  profile.  Also,  differences  in  bandwidth  performance 
between  the  types  of  delay  line  have  been  resolved  and  the  "buried- 
antenna"  design  has  been  shown  to  be  superior.  Ray  tracing  results  have 
been  correlated  with  an  experimental  result  from  each  type  of  delay  line. 

Delay  line  performance  is  typically  characterized  by  a  time- 
bandwidth  product  which  delineates  the  trade-off  between  bandwidth 
and  total  delay  time.  Magnetoelastic  delay  lines,  however,  should  be 
capable  of  large  time-bandwidth  products  and,  realistically,  are  not  well 
characterized  by  this  quantity.  The  bandwidth  of  magnetoelastic  delay 
lines  is  apparently  relatively  insensitive  to  the  frequency  range  of 
operation  and,  since  delay  times  are  directly  proportional  to  frequency, 
design  of  delay  lines  with  large  time-bandwidth  products  should  be 
possible.  This  should  be  especially  true  at  high  frequencies  and  for 
delay  functions  with  low  values  of  the  dispersive  delay  rate,  D. 

Consequently,  linearly  dispersive  magnetoelastic  delay  lines  can  now 
be  designed  with  an  a  priori  knowledge  of  the  attainable  bandwidth  and 
time-bandwidth  product,  and  the  required  magnetic  field  profile.  As 
described  in  Chapter  V,  bandwidth  is  maximized  for  magnetic  field 
profiles  which  cause  the  focusing  paramett.r  Q  t''  be  restricted  to  a 
limited  range.  For  a  given  dispersive  delay  function,  the  value  of  Q 
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is  determined  by  the  magnetic  field  profile  nonlinearity  constant  C,  and 
a  closed  form  expression  for  C  as  a  function  of  Q  is  available.  The 
selection  of  C  in  turn  allows  direct  estimation  of  attainable  bandwidth 
and,  ultimately,  time-bandwidth  product. 

Validation  of  the  ray  tracing  approach  comes  from  previous 
experimental  work.  In  1983  Borgeaud  [3]  constructed  a  "buried-antenna" 
delay  line  prior  to  the  development  of  the  described  ray  tracing  tools. 
His  selection  of  a  nominal  field  profile  constant  was  motivated  by  the 
simple  requirement  to  demonstrate  that  nonlinear  magnetic  fields  could 
produce  linear  dispersion.  Subsequent  bandwidth  analysis  of  his 
experimental  work  (a)  is  validated  for  his  experimental  case  and  (b) 
shows  that  his  choice  of  field  profile  constant  was  slightly  off  optimum. 
Earlier,  Platzker  [2]  had  built  and  tested  a  Strauss  type  delay  line  with 
synthesized  magnetic  field.  He  selected  a  linear  field  profile  (C  =  0), 
also  without  any  a  priori  method  of  estimating  bandwidth.  His  measured 
bandwidth  matches  ray  tracing  analysis,  and  a  larger  bandwidth  could  have 
been  attained  with  a  different  value  of  C. 

The  differences  in  bandwidth  performance  between  the  "buried- 
antenna"  and  Strauss  type  delay  lines  can  also  be  explained  by  ray 
tracing.  The  Strauss  type  delay  line  operates  over  a  narrower  band 
because  low  end  frequencies  cannot  reach  crossover  points.  This  effect 
is  apparently  the  sole  factor  in  Strauss  type  lines  exhibiting  reduced 
bandwidths . 

The  continued  analysis  of  magnetoelastic  delay  lines  requires 
further  correlation  of  ray  tracing  with  experimental  results.  Although 
the  validity  of  ray  tracing  in  magnetoelastic  media  is  well  established, 
predictive  properties  should  be  based  upon  a  variety  of  experimental 


cases.  Cases  of  special  interest  include  very  low  values  of  D,  the 


dispersive  delay  rate,  for  which  extremely  large  bandwidths  are  expected; 
values  of  C  close  to  zero;  and  values  of  C  very  large  negative. 
Additionally,  the  bandwidth  dependence  on  frequency  range  should  be 
exa.tiined  for  X-band  and  above.  The  case  of  D  very  small  is  of  special 
interest  since  previous  work  [31]  has  shown  that  a  constant  delay  profile 
is  when  focusing.  As  discussed  in  Section  III.E,  a  limit  on  D  for  Q  >  0 
car.  be  defined;  it  therefore  should  be  possible  to  design  delay  lines 
with,  relatively  small  fractional  changes  in  dispersive  delay  time. 

A  desirable  upgrade  to  the  existing  capability  would  be  the  ability 
;  r<=dict  attenuation  in  the  delay  ine .  The  existing  computer  program 
'r.rains  an  attenuation  calculation  subroutine  (adapted  from  Addison 
t:T  ai.  [23]).  There  are,  however,  very  few  experimental  data  to  confirm 
quality  of  the  approximations. 

.'■‘lysical  improvements  to  the  "bur  ied-a.ni  cnna"  magnetoej  ast  ic  delay 
i;r.e  ccnstructicn  are  also  desirable.  These  improvements  would  provide 
some  term  of  shielding  between  the  antenna  wires  to  improve  isolation, 
/'.itc,  adding  a  conducting  sheath  around  the  YIG  crystal  between  the 
r.  t  e  n a  holes  would  provide  a  "waveguide"  in  which  electromagnetic 
radiation  would  experience  cutoff.  Antenna-antenna  isolation  would  be 
.reproved  with  minimal  effect  or.  magnetoelastic  wave  propagation. 


Appendix  A 


MAGNETOELASTIC  WAVE  THEORY 


A . 1  Complete  Magnetic  Spin  Wave  Theory  and  Dispersion  Relations 

Magnetic  spin  waves  in  f  erriinagnet  ics  are  the  result  of  the 
magnet isa^ ion  field  and  its  incorporation  into  the  constitutive 
relations.  This  magnetization  field  is  induced  by  the  precession  of 
electron  spin  moments  about  an  applied  axis  of  magnetization.  When  a 
strong  enough  magnetic  field  is  applied,  all  the  magnetic  moments  precess 
around  .same  axis  of  static  .magnet  izat  ion . 

M.=.xv.’'-- 1  i  '  s  Equations  (in  MriS  units)  for  a  f  err  imagnet  ic  are  those  of 
a  sourceiess,  materially  isotropic  medium: 


V‘F  -  0 

(A. 

.1) 

V»F  =  0 

(A. 

.2) 

V  »  H  =  d  5  /  t 

(A. 

3) 

7  ■  ii  =  -aE/J^t 

(A. 

A) 

Tiie  constitutive  relations  vary  frotti  the  isotropic  case  in  that  the 
precr-sf  icn  of  magnetic  mv;m<^r,t'j  produces  time  varying  magnetic  fields 
each  of  which  generates  a  torque  on  its  respective  magnetization  vector. 
Thj.s  introduces  a  time  dependence  into  the  constitutive  relations 
(Figure  2.2); 

E  =  p, (H  +  M)  (A. 5) 

CM/ at  =  /p,/M  X 


(A. 6) 
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W.  =  cubic  anisotropy  energy 


_L  .  !fi 

(“o  '■)“  ^’'k 

magnetoelastic  coupling  constant 


elastic  displacement 


unit  vector 


r.cua;  ;ci 


Nd*M  =  demagnetizing  tensor  (geometry  dependent) 

a  -  9M  ,  ,  -  ,  j  • 

=-  -p^  M  »  ~  a  =  phenomenological  damping  constant 


A .  0  defines  the  interaction  between  the  internal  magnetic  field 


r.-u'r.'^;  izaticn.  H 


the  effective  internal  magnetic  field 


■.■k'r.’ a  1  terms  which  act  to  produce  torques  in  the  magnetization 


•an  be 


modeled  as  effective  magnetic  fields: 


--  the  applied,  external  static  magnetic  field; 

-  the  time  varying  component  of  applied  magnetic  field; 

--  a  phenomenological  damping  terra  which  accounts  for  losses. 

..  .  i  j  ■ 

"h.;:  'err:!  is  ignored  in  the  further  derivations  as  it  can  be 
!,r.;  eh  for  with  a  straightforward  addition  of  an  imaginary 
wave-.'ertcr  component  in  the  final  dispersion  relations; 

the  quantum  mechanical  exchange  force  has  no  analog  in 
f  ia:  .':;;ai  physics.  This  force  acts  to  align  nearest  neighbor  spin 
rrioments  and  is  required  for  ferro-,  antiferro-,  and  f  err  imagnet  ic 
t  f  ■  i  1  a  V  ;  - .  r  ; 
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(e)  --  demagnetizing  magnetic  field  caused  by  the  occurrence  of 
volume  divergence  and  surface  dipole  magnetic  terms  in  finite 
ferrite  samples; 

(f)  --  the  m.agnetic  field  induced  by  the  elastic  strain  of  the 
lattice.  This  term  can  be  treated  as  negligible,  although  it  is 
the  mechanism  for  coupling  between  spin  and  acoustic  waves.  The 
coupling  will  bo  introduced  subsequently  through  a  coupling  of  modes 
metliod ; 

(g)  --  the  anisotropy  field  which  depends  upon  the  orientation  of  the 
external  magnetic  field  with  respect  to  the  crystal  lattice 

St  ructure . 

Morgenthaler  has  developed  the  full  thecry  incorporating  all  these 
magnet oe 1  as* ’ c  terms  [20]. 

.h’Gglec;  i.ng  i)';e  smaller  terms  .'eav''.s  t.he  si.mplified  constitutive 
r  t;  i  a  t  i  o  r. : 


(  M  X  (H 


ext 
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(  A .  8  ) 


Assum.ir.i;  that  (a)  the  DC  cor.p'r.ent  r-  of  the  magnetization  and 
m.icr.ot'^-  fields  are  z-directed,  and  (b)  the  DC  components  are  much  larger 
than  tl'.o  time  varying  components; 


ft  -  M  z  +  m  |k!  =  M. 

i;  =  ii  z  +  h.  Ifli  =  I!,  t;  =  ti  ,  +  A  V^m 

r  f  ex 


( A  .  9  ) 


then  substitution  of  the  above  into  the  constitutive  relation  yields  the 
small  signal  constitutive  relation  which  explicitly  defines  the  relation 
between  the  time  varying  components  of  M  and  R: 


—  =  X  -  c^^m  -  m) 


(A.  10 
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Assuming  a  slowly  varying  plane  wave  nature  for  the  time  varying 
components : 


((-.t-k*r) 

^j(c>t-k*f) 

G 

(ot-k‘r) 


j  ( (ot  -  k*  r  ) 


and  equating  the  differential  operators  with  their  equivalents: 


V  -  i  k 


(A.  11 


(A.  12 


then  the  tensor  relation  between  the  small  signal  magnetization  and  the 
local  magnetic  field  can  be  defined.  The  coefficient  matrix  x  is  the 
Polder  Tensor. 


m  =  x*h 


(A. 13: 


The  small  signal  magnetic  flux  is: 


b  =  lip(h  +  m)  =  (A. 14 

M  =  MpCx  +  I) 

Since  the  relation  defines  the  extent  of  an  induced  magnetization  as  a 
function  of  applied  field,  the  Polder  Tensor  terms  are  susceptibilities. 

Applying  the  plane  wave  approximation  to  Maxwell's  Equations  in  the 
ferrimagnet ic  yields: 


-jk*te  =  0 

(A.  15 

-jk»b  =  0 

(A.  16 

-jk  >  h  =  jcjcg 

(A.  17 

-jk  *  e  =  -jwb 

(A.  *8 

ihe  F\F  magnetic  and  electric  fields  can  be  defined  as  a  function  of  the 
small  signal  magnetization  field: 

-k(k*m)  +  kpfn 
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(A. 20) 


and  substitution  of  A. 19  into  A. 10  yields  a  complete  magnetic  spin  wave 
dispersion  relation  (Figure  A.l): 


k  = 


(p^-p-x^)sin^(e)+2p  ±  {  ) sin^  (  6 )-t-4x^cos  ^ ( 9 ) } 
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The  dispersion  relation  provides  for  an  ordinary  wave  with  "EM-like" 
characteristics  and  an  extraordinary  wave  which  is  the  magnetic  spin 
wave.  The  extraordinary  wave  is  right-hand  elliptically  polarized. 


A . 2  The  Magnetostatic  Approximation 

From  Equations  A.  19  and  A. 20,  for  the  case  |kj  >>  k^^,  the  magnetic 
and  electric  fields  are  approximately; 

-k(k*ffi)  ■*  k^m 
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The  electric  field  is  proportional  to  1/k  and,  as  (k|  becomes  very  large, 
a  magnetostatic  approximation  can  be  made: 


e  ^  0 


(A. 24: 


V  X  h  =  0 


(A. 255 


The  magnetic  field  can  be  modeled  as  the  gradient  of  a  scalar  magnetic 
potential  field  and  a  simplified  approximation  to  hj.£  thus  becomes 
available : 


(A. 26) 


(A. 27) 


(A. 28) 


Substituting  Equation  A. 28  into  Equation  A. 10  produces  an  equation  which 
is  homogenous  in  in  and  whose  secular  determinant  gives  the  dispersion 
relation  for  an  infinte  ferrite  under  the  magnetostatic  approximation; 


=  Wj,jSin^(G)) 


(A. 29) 


The  existence  and  propagation  of  elastic  waves  in  isotropic  YIQ  are 


Having  separately  derived  the  equations  for  the  magnetic  spin  wave 
and  elastic  wave  systems,  the  coupling  of  modes  approach  can  be  used  to 


define  their  interaction.  This  treatment  is  equivalent  to  including  an 


hjjjg  term  in  Equation  A.  5.  The  complete  approach  has  been  extensively 
developed  and  yields  the  approximate  dispersion  relation  [7]: 
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Conversion  between  spin  and  elastic  waves  occurs  when: 


(A. 33) 


cj  =  -  Vg^  =  elastic  wave  velocity  (A.  34) 

^ew 

Although  the  splitting  term  is  proportional  to  frequency,  it  is 
generally  of  an  order  of  lO’  Hz  at  room  temperature,  typically  less  than 
1%  of  the  unperturbed  frequency  value.  The  dispersion  relation  is  shown 
in  Figures  2.8  and  2,9. 

The  efficiency  of  coupling  is  determined  by  the  ratio  of  local 
magnetic  field  slope  to  a  critical  value  [20]: 
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(A. 35) 


As  pointed  out  by  Wadsworth  [33],  this  condition  is  easily  met  and  is  not 
a  good  design  criterion. 
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Appendix  B 


MAGNETIC  FIELD  REQUIRED  FOR  LINEAR  DISPERSION 


The  tirae  spent  by  magnetoelastic  waves  propagating  through  the  delay 
line  can  be  accounted  for  by  two  wave  behaviors:  (a)  spin  exchange  wave 
and  (b)  elastic  wave.  Transitions  between  wave  types  are  assumed  to 
occupy  negligible  time.  The  backward  magnetostatic  wave  contribution  to 
total  propagation  time  is  also  neglected  because  the  wave  has  very  high 
group  velocity. 


TD(f)  =  T  +  T 

^  sew  ew 


The  tirae  spent  as  an  elastic  wave  which  has  fixed  velocity  is: 


(B.l) 


(B.2) 


The  group  velocity  of  the  spin  exchange  wave  is  a  function  of  position 
since  the  raagnetic  field  is  inhomogenous .  The  time  spent  as  a  spin 
exchange  wave  is  calculated  as  follows: 


-n  „  cp  dz 

^sew  ~  ^  V 


(B.3) 


Differentiating  the  approximate  spin  exchange  wave  dispersion  relation  of 
Equation  2.22  gives: 


(B.4) 


and,  for  a  spin  exchange  wave,  the  axial  group  velocity  is 
(Equation  2.23): 


z  M  ex  z 
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Substituting  Equations  B.4  and  B.5  into  Equation  B.3  gives  a  solvable 
integral  in  the  variable  instead  of  z: 


T  =2 
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l^cp^  l^tpl  ^“d,  from  Equation  2.28,  the  spin  wave  propagation 
time  can  be  estimated: 
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sew 


H'(Zcp) 


(B 


The  magnetic  field  at  the  crossover  point  can  be  explicitly 
calculated  from  Equations  2.22  and  2.27: 


'^0  ew 

The  second  terra  represents  a  very  sraali  contribution  to  the  field 
strength  at  crossover  (-1%)  and  can  be  neglected  (as  done  in 
Equation  3.5).  The  field  strength  at  crossover  can  therefore  be 


considered  to  be  proportional  to  frequency  and  the  total  time  spent  is 
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From  Equation  B.9,  the  spatial  derivative  of  ^^cp^  thought 


of  as  proportional  to  a  spatial  derivative  of  f,  and  substituting  into 
Equation  2.10  with  Equations  B.l  and  B.2  gives  an  exact  differential: 
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(B.12) 


Equating  this  function  of  frequency  to  the  desired  linear  delay 
function  in  Equation  3.1  and  integrating  gives: 
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It  is  apparent  that  appropriately  locating  the  crossover  point  z 


cp 


as  a  function  of  frequency  allows  the  linear  dispersion  to  be  obtained. 
The  field  strength  at  crossover  is  directly  proportional  to  frequency 
(Equation  B.9)  and  a  straightforward  substitution  into  Equation  B.15  can 
be  made: 
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Rearranging  Equation  B.16  and  solving  the  quadratic  for  H(z=z  )  yields 
an  expression  for  the  magnetic  field  strength  required  for  the  linear 
delay : 
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Generalizing 
derivation  of  the 


z^p  to  cover  all  values  along  the  z-axis  completes  the 
required  one-dimensional  magnetic  field. 
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The  fixed  delay,  T^,,  must  be  negative  for  the  magnetic  field  to  be 
positive  everywhere  in  the  rod.  A  change  of  sign  in  the  magnetic  field 
implies  that  domains  might  form. 

The  magnetic  field  prescribed  above  is  for  z  >  0,  where  z  =  0  is 
the  midpoint  of  the  rod.  Replacing  z  with  |z|  gives  a  function  which 
is  symmetric  in  z  but  has  discontinuous  derivatives  at  z  =  0 .  Since 
continuity  in  H(r,z)  and  its  derivatives  inside  the  rod  is  required 


by  Maxwell's  Equations,  an  approximation  near  z  “  0  is  required.  A 


Appendix  C 


MAGNETIC  FIELD  FOR  A  YIG  ROD 


Maenetic  Field  in  a  YIG  Rod 


Ray  tracing  in  the  magnetoelastic  delay  line  requires  knowledge  of 
the  static  magnetic  field  everywhere  within  the  YIG  rod.  Morgenthaler 
has  devised  a  method  for  estimating  this  magnetic  field  by  expanding  the 
specified  on-axis  field  subject  to  certain  constraints  [1]. 

Maxwell's  equations  for  static  fields  inside  a  f err imagnet ic  yield 
the  following  constraints: 


M  x  H  =  0 


M*H  >  0 


M  =  M, 


V  X  H  =  0 


V»B  =  0 


Since  V  x  H  is  zero,  both  H  and  M  can  be  expressed  as  the  gradient  of  a 
scalar  field: 


H  =  -Vi// 


M  =  M 

ivV;|  ° 


The  scalar  field  can  be  expanded  as  a  series; 


Although  H  and  M  could  be  expanded  in  terms  of  the  same  coefficients, 
they  can  also  be  expressed  as  follows: 


00  00 

2n  .  =  V  n/_  .  /_\_2n+l 

2n' 

n=0  n=0 


H(r,z)  =  z  y  a,^(z)r'^'^  +  t  ^  2(n  +  l)32n+l^^^ 


(C.9) 


M(r,z)  =  M, 


2  I  +  t  I  b,„^,  (z)r^”'^^ 


n=0 


n=0 


(C.IO) 


The  coefficients  a„  and  b„  are  functions  of  z  and  can  be  determined  by 
n  n  •' 

the  following  expressions: 


b  j  (  z ) 


I  " 

-j 

a] (z)  +  Mg 


(C.ll) 


a^Cz) 


b^  (z) 


2  bj(z)aj,(z) 


-i  b^z) 


(C.12) 

(C.13) 


-1  ,  " 


(39<=‘*'^nb7o)  +  1  (2(s-k+l) 


2(s+1)^“2s"“o'^2s 


k=l 


^2(s-k+l)^2k  ^2k‘^2(s-k)  +  l  ^ 


'2s  +  l 


Sn  + 

0  0 


(C.14) 


‘2(s+l) 


®o^2s  +  l  ^  t^^®'‘^'^^^®2(s-k+l)‘^2k  ®2k^2 ( s-k)  +  l  ^ 
_ k-=l _ 

2(s  +  1) 


(C.15) 
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‘’2(3  +  1) 


^  \'^2s-k+2  2  '’s+1 

t  *• 


(C.16) 


The  solution  is  implemented  by  invoking: 


^  ‘ 

ajj(z)  =  H(z) 
a^Cz)  =  H' (z) 


d"(a(,(z). 

^  d”"‘H(z) 

dz" 

A 

dz 

a^Cz)  = 

r 

H(  z)dz 

(C.17) 

(C.18) 


(C.19) 

(C.  20) 


C . 2  External  Potential  Field 

C.2.A  Outer  Field  Expansion 

Morgenthaler ' s  method  for  synthesizing  a  desired  magnetic  field  in  a 
rod  or  disc  requires  the  calculation  of  equipotential  surfaces  outside 
the  rod  [1].  Once  the  magnetic  field  expansion  for  the  interior  of  the 
rod  is  known,  an  expression  for  a  scalar  potential  field  outside  the  rod 
can  also  be  found.  From  this  scalar  potential,  the  required 
equipotentials  can  be  found. 

The  magnetic  scalar  potential  field  outside  the  rod  must  be  a 


solution  to  the  Laplacian  in  cylindrical  coordinates: 


The  conventional  selection  of  Bessel  functions,  however,  becomes 


difficult  to  solve.  Alternatively,  the  potential  field  can  be  expanded: 


=  I  +  I  c'Q*(r,z) 


ext  ^  n  n 
n=0 


(C.22) 


and,  with  solutions  using  a^Cz)  =  z  ,  a  linear  sum  of  mixed  polynomials 
each,  of  which  is  a  solution  to  the  Laplacian  can  be  derived: 


(C.23) 


3  2 
2 


Sz'r^  +  I  r' 


rhese  polynomials  resemble  the  Legendre  polynomials.  A  second  set  of 


1  yr.om  1  a  1  s  is  also  included: 


=  ln(r)  +  1 


(C.24) 


z  ln(r)  +  1 


ln(r)  ~  2  in(r)  +  z^ 


z^  ln(r)  -  —  zr^  ln(r)  +  z' 
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Q  =  ln(r)  -  3z^r^  ln(r)  +  |-  ln(r)  +  z* 
’  o 


Recursion  relations  which  explicitly  define  these  functions  are  described 
by  Morgenthaler  [1]. 

* 

C . 2 . B  Boundary  Conditions 

The  coefficients  in  the  outer  field  expansion  are  found  by  solving 
the  boundary  conditions  at  the  ’'od's  cylindrical  surface.  The  boundary 
conditions  for  magnetic  fields  are  continuity  of  (a)  tangential  H  and 
(b)  normal  B.  The  outer  field  potential  can  then  be  expressed  cirectly 
as  a  function  of  the  inner  field. 

In  order  to  simplify  the  calculation,  the  internal  tangential  H  and 
normal  B  fields  are  approximated  by  polynomials  in  z  at  the  surface.  The 
conversion  from  the  nonlinear  expansion  inside  the  rod  to  a  polynomial  at 
the  inner  surface  of  the  rod  uses  an  integral  mean  square  error 
minimizat ion : 


E  = 


n=0 

length/2  ( 
-length/2 


(C.25) 


n«0 


n=0 
2 


dz 


length/2 

N 

bJ'  '(R,z) 

-  1 

-length/2 

n=0  ^ 

dz 


(C.26) 


Bj^'^^(R,z)  =  hJ"''(R,z)  +  Mi''’^(R,z) 
R  =  rod  radius 
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Minimizing  the  error  with  respect  to  each  set  of  the  polynomial 
coefficients  gives  N  linear  equations  of  the  form: 


length/2 

-length/2 


2  hJ’'^(R,z)(^)^ 


(-)"  I  «  i-T 

n=0 


(C.27) 


The  first  term  in  each  integral  is  solved  by  normalizing  the  length  to 


[-1,1]  and  applying  a  Gauss-1 egendre  series  approximation: 


1  N 

f(u)du  -  Y.  f(.u^)w. 
-1  i  =  0 


(C.28) 


u.  =  zeroes  of  Gauss-Legendre  function  of  order  N 
^  weights  of  Gauss-Legendre  function  of  order  N 


The  second  term  is  directly  integrable; 


length/2  N  N  a  L  .  , 1  length/2 

Y  ^  rz-\n+i  _  Y  ^  /■z\n+i+]  | 

,  ^  ^  ^  n  +  i  +  1  -length/2 

-length/2  n=0  n=0 


(C.29) 


length/2  N  N  /5  L  .  ,  length/2 
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■length/2 


2/JnL  ,  .  , 

'  ^  /■  1  \n-t-x+l 


=  Y  .  n  +  i  +  1  '•2^ 


n-fi  +  1  odd 


n+i+1  even 


I 


The  coefficients  and  fi ^  are  found  by  equating  the  two  expressions  and 
solving  two  systems  of  N  linear  equations.  This  approximation  to  the 
internal  field  at  the  surface  is  readily  found  by  computer. 

Calculation  of  the  external  field  requires  fulfilling  the  boundary 
conditions  at  the  rod  radius: 

«  H^"’^(R,z)  -  h|^*^(R,z)  (C.30) 

=  b|^^(R,z)  (C.31) 

The  mean  square  error  at  the  boundary  can  be  expressed  as: 
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N 

-  I  C.f. 

j  =  l 


1  + 


I  Vii- 


int 


30 


int 


3f. 


^  2 


3r 


J  =  1 
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(C.32) 


With  the  inner  field  defined  as  a  polynomial  in  z,  the  least  squares 
method  can  be  applied.  An  error  expression  is  formed  and  differentiated 
with  respect  to  each  of  the  outer  field  coefficients  (with  equal 
weighting  of  the  boundary  conditions  in  the  error  expression), 
.'iorgenthaler  and  Platzker  demonstrated  the  closed  form  solution  for 


the  outer  field  coefficients  [2]: 


(C.33) 


^  «  V2k.k^Cp-2k  "  ^  ^k)S^2kJ 
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Rft^  N-p 


C  = 
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S  --=  ° 


^.k  = 


n !  4 


«k  = 


(n  -  2k) !  (k! )■ 

k  =  0 


1  -  I  7 


k  >  0 


i  =  l 


(C.35) 

(C.36) 


f 

The  procedure  involves  assuming  that  Cj^  =  0,  and  subsequently  finding 

t  • 

C  i,C  o,C  T,  ...,C.  Once  these  coefficients  have  been  found,  the 
external  potential  field  and  its  gradient,  the  magnetic  field,  are 
de f  i ned . 

Wadsworth  has  noted  an  interesting  correlation  between  external 
field  and  delay  function  [33],  He  found  that  the  coefficients  for  delay 
function  1000  f  -  200  were  approximately  twice  those  of  the  delay 
function  500  f  -  200.  Doubling  the  external  field  coefficients 
corresponds  roughly  to  doubling  the  energizing  field. 
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